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ABSTRACT 

In this paper the special theory of relativity is extended to include the vectorial and scalar 
transformations of kinematical and dynamical quantities associated with bodies whose proper or 
rest mass is variable with time (i. e. rockets). Besides the known transformation formulas for 
force (or momentum flow rate) the general transformation formulas for power (or total energy flow 
rate) are given for the first time. All text books on the theory of relativity introduce the classi- 
cal definition for the power Which is an overspecification. It would be correct only when the rest 
mass is not changeable with time. This assumption is always fulfilled when applying the theory 
of special relativity to fast moving electrons, atoms, or nuclear particles, but it is not true for 
fast moving rockets. Generalized relativistic conservation laws of momentum and total energy 
(mass) are derived. 

An application of relativistic dynamics to rocket propulsion gives the data of an arbitrarily 
accelerated rocket in free space (without external forces) in the system of a stationary earth 
observer, and in the rest system of an astronaut centered in the moving rocket itself. 

Two special cases of rectilinear motion of a rocket with constant exhaust velocity are 
treated: 


(1) constant thrust acceleration (hyperbolic motion) 

(2) constant mass flow rate or constant thrust 

Tables with numerical values for dimensionless flight parameters will be given for both cases. 


In this report Einstein’s general theory of relativity (gravitational theory) is applied to the 
motion of an artificial satellite revolving in an arbitrary orbit around a central body and the time 
dilatation effect for this satellite is given. This relativistic perturbation theory is based on 
Einstein’s general field theory, differential geometry of non-Euclidean spaces, potential theory, 
and celestial mechanics. The short periodic perturbations are excluded by using time average 
values over a revolution. The secular and long-periodic (non-relativistic as well as relativistic) 
perturbations of the osculating orbital elements, which represent deviations from the elliptic orbit, 
are presented here for the case of a rotating, non -homogeneous, oblated spheroidal central body. 
This is an extension of the work of Einstein (1915) who considered motion around a mass point 
as well as the work of deSitter (1916) and, independently, of Lense and Thirripg (1918), who 
treated the relativistic motion around a rotating, homogeneous, spherical central body, omitting 
the terms due to the square of the angular velocity. 

A formula for the relative difference of the time rates of a satellite clock, compared against 
a standard earth clock (time dilatation effect) is derived for orbits of any eccentricity and equa- 
torial inclinations, thus extending the paper of Winterberg (1955), Singer (1956) and Hoffmann 
(1957). 
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ASTRO REL ATI VI TY 


By 

Helmut G. L. Krause 


PART I 

RELATIVISTIC ROCKET MECHANICS 


SUMMARY 

In this part the special theory of relativity is extended to include the vectorial and scalar 
transformations of kinematical and dynamical quantities associated with bodies whose proper or 
rest mass is variable with time (i.e. rockets). Besides the known transformation formulas for 
force (or momentum flow rate) the general transformation formulas for power (or total energy flow 
rate) are given for the first time. All text books on the theory of relativity introduce the classi- 
cal definition for the power which is an overspecification. It would be correct only when the rest 
mass is not changeable with time. This assumption is always fulfilled when applying the theory 
of special relativity to fast moving electrons, atoms, or nuclear particles, but it is not true for 
fast moving rockets. Generalized relativistic conservation laws of momentum and total energy 
(mass) are derived. 

An application of relativistic dynamics to rocket propulsion gives the data of an arbitrarily 
accelerated rocket in free space (without external forces) in the system of a stationary earth 
observer, and in the rest system of an astronaut centered in the moving rocket itself. 

Two special cases of rectilinear motion of a rocket with constant exhaust velocity are 
treated: 

(1) constant thrust acceleration (hyperbolic motion) 

(2) constant mass flow rate or constant thrust 

Tables with numerical values for dimensionless flight parameters will be given for both cases. 



RELATIVISTIC ROCKET KINEMATICS 


A. General Lorentz Transformations 

The special theory of relativity is based on two principles: 

1. The postulate of relativity. It is impossible to measure or detect unaccelerated translatory 
motion of a system through free space. 

2. The postulate of the constancy of the velocity of light. The velocity of light in free space 
is the same for all observers, independent of the relative velocity of the source of light and the 
observer. 

Using these postulates, A. Einstein (Ref. 1) derived the Lorentz transformations (trans- 
formation formulae between two reference frames in relative uniform motion to each other) previ- 
ously derived by H. A. Lorentz (Ref. 2) and H. Poincare (Ref. 3) on the basis of the electron 
theory. 

Two coordinate systems S and S' moved against each other with a constant system velocity, 
(where v is the velocity vector of the system S ' relative to S and v ' = - v is the velocity vector 
of the system S relative to SO have origins 0 and O', respectively, coinciding at the time t=t =0. 
The transformation of a position vector r = (x, y, z) and the time t in the system S to the corres- 
ponding quantities r' = (x ', y z') and t'in the system S ' is given by the general Lorentz trans- 
formation (without rotation) in a paper by G. Herglotz (Ref. 4), namely 



where* c is the light velocity and 

a * >J 1 - ( v/c) 2 <1 ; 1 — a 2 = ( v/c) 2 (3) 

All these transformations satisfy the equation 

s 2 = (r • r) — c 2 t 2 ~ (r'« x') - c 2 t' 2 ( invariant) (4) 

The inverse equations (2) follow from equations (1) by interchanging (r' , t ) and (r, t) 
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and replacing v by v — V. 

The general Lorentz transformation with rotation (when the Cartesian axes in S and S' do 
not have the same orientation) are 



where D is the rotation operator (a tensor or matrix). Due to 

0v' = - v ; v'= - D 1 v 


there is 



and the inverse relations 



because of the identity 



(5) 


( 6 ) 

(7) 


( 8 ) 


(9) 


These equations also satisfy equation (4) . For inhomogeneous Lorentz transformations (where 
the origins 0 and 0 ' do not coincide at the time / = t ' = 0 ) the distance s in the space-time world 
will no longer be invariant. The invariant quantity for these general transformations is now the 
square of the four-dimensional line-element ds 2 = dx 2 + dy 2 + dz 2 - c 2 dt 2 . 


Special Lorentz transformations (where the constant velocity v is in the direction of the 
positive x- axis, thus v = [ v, o, oj) are 


x ' = — = x cosh \fj — c t sinh \jj 


y = y ; z = Z 
c t — x v/c 


c t = 
where 


- c t cosh — x sinh xjj 


X + Vt , . . , . . , 

* = = x cosh xjj + c t sinh xjj 


y = y , z - z 

. ct'+ x' v/c 
C t — — 


- c t' cosh if/ + x' sinh if/ 


if/ = cosh^ 1 (1/a) = cosh _1 — = tanh -1 (v/c) 

vr- (v/c)* 


( 10 ) 


3 



B. Transformation of ParticU Velocities 


The motion of an arbitrarily moving particle will be given in the system S by its instantane- 
ous position x = x (t), y = y (t), z = z (t) or r (t). The particle velocity vector is u « dx/dt =* 
* (u x ,u y ,u x ) = (£,y,£) and the velocity itself is given by u « (u 2 + u 2 + u 2 )L Corresponding 
quantities in the System S 'are marked by a prime . 

By differentiation of eqs. (1) and (2) one obtains the transformation formulae for the 
time elements and for the Lorentz contraction factor a: 



Transformation formulae also follow for the local velocity vector of a particle in the reference 
system S and in the system S' in relative uniform motion with constant velocity v to S, namely 

, dr' dr'/dt _ au4-v i (1 - a) (u»v) /v 2 ~ 1 j , 

U “ dt' " dt'/dt l-(u*v)/c 2 


dt dx/dt* a u ' 4- v { (1 - a) (u '»v)/^ 2 + I 1 
dt = dt/dt' ~ 1 + (u '*v)/c 2 


(14) 


These formulae are the basis for the relativistic kinematics. By squaring eqs. (13) and 
(14) there results 



(15) 

(— ) 2 - 1 - a 2 1 ~ (u 1 ; u4c 

(16) 

combining eqs. (jl)> (15) and (16) yields 

dt' 1 — (u*v)/c 2 a ll-(u/c) 2 

dt a - l + (uW)/c2 - V l-(uVc)* 

(17) 

and thus 

dtyj\—(u/c f" = dt' \j\ ~ (u'/c)*' = dr (invariant) 

( 18) 


where r is called the proper time measured in a system centered in the moving particle. 


For u j| v respectively u'jjv eqs. (13) and (14) yield the usually quoted addition theo- 
rem of particle velocities 


u 


U - V 


1 - («J •v)/c2 


u '+ V 

° = l + (u '.y)/c2 


(19) 
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For uj_ v, respectively, ul v eqs. (13) and (14) yield 

u '= au - v ; u = au '+ v 


( 20 ) 


For [u| = c (photons) there is u ' = - u, while for u = v there follows u ' = 0 (transformation to rest) 
meaning that the coordinate system S ' is centered in the particle itself (rest system S ). 

For the special case where v is in the direction of the positive x - axis, eqs. (13) and 
(14) give 


l-(u x v/cV ' y l-(u x v/c 2 ) ' 


a u y 

Uy 1 +(u'v/c 2 ) 


x 1 + (u'v/ c 2 ) 

Taking scalar and vector products of eq. (13) with v produces 
, u*v — V 2 


au z 

(21) 

1 — (u x v/ c 2 ) 

aU z 

1 +(v x v/c 2 ) 

(22) 


1 — (u*v)/c 2 


a [u x v| 

1 — (u*v)/c : 


(23) 


Introducing the angles 0 = (u, v) and Q' ~ < (u \ v) the above equations can be written after 

division by v, in the following form 


U ' COS 0 ' — 


u cos 0 — 


1 - 


U V COS 0 


u ' sin 0 ' = 


a u sin 0 


1 - 


U V COS 0 
C2 


(24) 


and thus 


[ (^ cos 6 — v) 2 + ( a u sin Q) 2 ] */2 
u v cos 0 

1 C 2 


and 


tan Q' = 


a sin 0 
cos 0 — ( v/u ) 



The inverse equations follow at once by interchanging the primed and unprimed quantities and re- 
placing v by — v y thus 


u cos 0 


u ' cos 0 ' + 


1 + 


u 'v Cos 0 ' 


u sin 0 = 


au ' sin 0 ' 


1 + 


u 'v cosQ ' 


(27) 


_ [(a / cos O' + v) 2 + Kg u ' si n Q Q 2 j 1 / 2 _ , [1 + 2 ( v/u ') cos Q ' + ( v/u 0 2 - ( v/c) 2 sin 2 Q '] V 2 
, u 'v cos O' _ u'v cos 6 

1 + o 1 + Ji 


(28) 


tan 0 = 


a sin 0 ' 
cos 0 ' + ( v/u ') 


(29) 
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Applying eqs. (24) and (26) respectively eqs. (27) and (29) to photons (u = u' - c) 
the relativistic formulae for the aberration of light follow: 


cos 0 -(v/c) . a sin 0 a sin 0 

cos Q = -5 — /. i ■-■■■ i a * ; 0 - t — ; tan <9 = 


1 - ( v/c) cos 0 


1 — (v/c) cos 0 


cos 0 — ( v/c) 


( 30 ) 


^ cos 0+ (v/c) . - a sin 6 ' n a sin 0 ' 

cosd= - ; «n*- ; tanfl = 


1 + ( v/c) cos 0 ' 


cos 0' + (v/c) 


(31) 


A more rigorous proof can be derived from the invariance of the phase of an electromagnetic 
wave. This principle also gives the relativistic formula for the Doppler effect. If v is the 
frequency, w the phase velocity, n the wave normal or the unit vector in the direction of the ray, 
and [k - ( v/w ) n] the wave propagation vector, the invariance of the phase can be expressed by 


v't'-(k'*r/) -vt — (k*r) or — v 


( 32 ) 


Eliminating r' and t' by means of eq. ( 2 ) and comparing space-dependent and time-dependent 
terms respectively on both sides of the above equation yields the transformation formulae 


k = k'+ v 



kW 

/ 

V 

LV« ) 

V* 

+ ac*j 


v'+ (k W) 


The elimination of r and t by means of eq. (11.2) leads to the inverse relations 
k' = k + v| 


\(L _ i) 

r 

* 

> 

, v-k‘v 

U y 

v 2 a c 2 

, v — 

c. 


( 33 ) 


( 34 ) 


Introducing the direction angle of the wave normal Q = < (n, v) the transformation formula for the 
frequency can be written 


, v — (k*v) 1 - (n»v) /w 1 - (v/w) cos 0 

v = — = v t — ~ v — — — — 


( 35 ) 


Taking the scalar and vector products of k ' with v produces 

, ✓ (k*v)-v (v/c) 2 v' i , \ v (n*v) - w (v/c) 2 

k »v = or — \n • v; = — 

a w w a 

and 

|k'x.v| = |k x v| or -^r|n'xvl = |n X v| 

These equations can be written 

v ' n , v cos 0 - (vw) /c 2 v' . v -a 

— ^ cos 6 = — — — — ; — > sint/ = — sin u 

w w a w w 


(36) 


Dividing both equations gives the transformation formula for the direction of the wave normal 

a sin0 


tan 6 ' = 


cos 0 - (v w)/c 2 


(37) 
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Squaring and then adding both equations produces 


— — [ sin 2 0 + {cos 0 — (vw)/c 2 } 2 /a 2 ] 1/2 
w w 

Eliminating v'/v yields the transformation formula for the phase velocity 

, 1 — ( v/w) cos 0 

W = W [a 2 sin 2 0 + {cos 0 — ( vw)/c * t *] 1/2 

w;— V cos 0 


1 f V • /5\ 2 o ^ 

1- (t y “ 2 cos0 + (tV 

The inverse transformation formulae for the wave characteristics are 


2 1/2 


( 38 ) 


, 1 + ( v/w ') cos 0 ' n a sin 0 ' 

v = v - ; tan 0 = 


cos 6' + (vw )/ c 2 


(39) 


, 1 + (v/ w ') cos 6 ' 

W W [a 2 sin 2 6' + {cos 0' + ( vw ')/c 2 1 2 ] 1/2 

w ' + v cos 0 ' 


i ( v • a '\ 2 o vw ' a' (vw'\ 

1 - \T 6 ) + 2 ~ cos 6 + V7 5 " ) 


2 1/2 


(40) 


A comparison of eq. (26) with eq. (37) and eq. (38) with eq. (25) shows that u = c 2 /w 
and u' = c 2 /«/. Therefore velocity and direction of a particle are transformed in the same way 
as corresponding quantities for a wave with phase velocity w - c 2 /u. DeBroglie used this re- 
sult in his wave theory of elementary particles. 


Applying the above mentioned formulae to a light source in vacuo (w — w' = c) the follow- 
ing transformation formulae result 


cos 0 ' = 


cos 0 — — 
c 

1 — — cos 0 
c 


1 / = V — 

V 


1 - ^ cos 6 


, 1 + cos 0' 


sin 0 ' 


a sin 0 

1 — — cos 0 
c: 


; tan 6 ' = 


a sin 0 


cos 0 - 


(41) 

(42) 


cos 6 ' + 


cos 0 = 


1 + — cos 0 ' 
c, 


sin 0 — 


a sin 0' 

1 -I- — cos 0 ' 
c 


; tan 6 — 


a sin 0' 


cos 0' + 


(43) 


tan 


0' /l + 


i-L 


0 t- 

tan _ ; 1 + - 


cos 0 ' — 


1 - — cos (9 
c 


(44) 


It is easy to derive the classical aberration formula from these equations, namely 
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sin {O' — 6) - sin 0 ' cos 0 - cos 0' sin 0 

- -(I -a) cos 0 

= sin 0 — 

1 - - cos 0 


(T - a) cos 0 ' + - 

= sin 0 ' c 

1 + - cos 0' Q'_ Q Sin 0 

~ r ' 

For — <5C 1 (a » 1) the above mentioned formulae yield 
c 

0'-0~ - sin 6' ~ -sin# , 
c c 

the same equation which follows directly from the figure. 

In discussing the Doppler effect there are two special cases 
a) the usual longitudinal Doppler effect: n || v {0 = 6'- 0) 



~ sin 0 


(45) 



a 

b) the transverse Doppler effect: n JL v 



77 - sin’ 1 a) 


(46) 


, v v 

a \j 1 - {v/~c)2 


> V 


(47) 


The decrease of the frequency, v ' , of the emitted light corresponds to a shift towards the red 
in accordance with the time dilatation of a moving clock. 


It is useful to have available the transformation formulae of Certain other quantities. Using 
eq. (44) the transformation formula for the solid angle d Q = sin 6 dO dcj) 


dQ' ^ d (cos O') 
d£l d (cos 6 ) 


a 


2 


(1 + - cos O') 
c 


2 


(1 - — cos 6) 2 

C 


a 


(48) 


The transformation formulae for the amplitude A, the volume V of a laterally bounded, finite 
wave and the total Energy E - V 2 A 2 V of the wave £re 


A ' 

7T 


v 

v~' 


E ' 

ir 


v 

V 


1 - - cos 0 
c 


1 + - cos 0' 
c 


( 49 ) 


The total energy density u = E/V = V 2 A 2 will be transformed according to 
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u 





c 


a 


cos 0 


2 



cm 

<m' 


thus 


( 50 ) 


u ' dQ, ' = u dQ 


Due to the transformation formula for the number density of photons D — u/hv : 


D' 


D 



(51) 


the number of photons N — D» V — E/hv is a relativistic invariant: N'= N. Thus the relativi- 

stic transformation from one to another coordinate system is not connected with a creation or 
destruction of photons. 


For practical applications of the relativistic Doppler effect and aberration formulae, it is 
more advantageous to replace the angle 0 between the velocity vector v of the light source 
(relative to the observer) and the actual direction from the light source to the observer by the 
angle 0 = 0 + n between v and the actual direction from the observer to the light source. Like- 
wise the angle 6 ' has now to be replaced by 0' = 6' + n for the corresponding apparent angles. 
Furthermore, v ' should be the proper frequency v 0 of a light source or a transmitter while v is 
the frequency measured by an observer or a receiver. The eqs. (41) to (44) therefore yield 
( /3 = v/c ; a = \/l - /3 2 ) 


1 +8 cos © 1-/3 cos 0 ' 

v = v , — ; v - v 5 


Vi -|S 5 


Vi -P 2 


cos © + B • yjl - 8 2 Sin® ^ , y/1 - sm© 

os © '= cL- ; sm © = A- L- ; tan© = ^ 

1 + 8 COS 0 1 + /3 cos 0 cos © + /3 


© = 


p 

cos Q /3 
1 — /3 cos © 7 1 


■p 


. ^ v/1 — B 2 sin ® ~ v - - 

m © = rS~v— ; tan © = — s - 

1 - /3 cos 0 cos 0-/3 


\f 1 - /3 2 sin © ' 


0 , 0 1 -B 2 

tan — = tan T ; 1 - ^ cos 0 = x ~ ff cos q 


0' /T^r 
2 “ Vi + B 


(52) 

(5*) 

(54) 

(55) 


0 = ©' = 0 (receding space vehicle): 


v/i v = — = 5- =/-! 77 < 1 (shift to the red) 

0 a 1 + £ "V 1 + P 


0 = -~- ; © ' = cos 1 jS: 


= \/l — f3 2 = a <1 (shift to the red) 


0 ' = ; 0 = 77 — cos ^ 

4 


1 ^/ 1 ^ = ^ = — > 1 (shift to the violet) 

* yl - p a 


0 = 0 ' = 77 (approaching space vehicle): 


v/v = l±JL = a 

a i - /S 


1 + B 

— > 1 (shift to the violet) 

1 ~ P 



The above mentioned formulae also apply to the case where the ( earth) observer is moving 
with respect to the fixed stars and the light source (fixed stars) is at rest. It is only necessary 
to interchange primed and unprimed quantities and to replace v by ~v. 

An illustration of the before mentioned formulae is shown in the following table: 


Table 1. Aberration and Doppler Effect 
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<r= c 

n -1 P+cosd 


v/v 0 

a- p 2 ) 1/2 


1 + jS cos 0 


1 + /3 cos 0 

0 

P = v/c 

0.5 

0.995 

0 

f3 = v/c 

0.5 

0.995 

0° 

0° 

0° 

0° 

1 

0.5773 

0.0501 

30 

30 

17.588 

1.537 

1 

0.6043 

0.0537 

60 

60 

36.869 

3-311 

1 

0.6928 

0.0668 

90 

90 

60.000 

5.732 

1 

0.8660 

0.1000 

120 

120 

90.000 

9.912 

1 

1.1547 

0.1990 

150 

150 

139.792 

21.166 

1 

1.5274 

0.7230 

180 

180 

180.000 

180.000 

1 

1.7321 

20.0000 


C. Lorentz Contraction and Time Dilatation 


Two important applications of the Lorentz transformations can be made: 

1. A measuring rod at rest (u « 0) relative to the system S has the length Ar«-r 2 — r x . In the 
system S ' moving with the velocity v relative to S the two end-points of the rod have a relative 
velocity u'»— -V and simultaneously 


[ A /'■* t' - t' = 0 or A / = (Ar*v)/c 2 » KAr*v)/t; 2 } (1 - a 2 )] 

have a distance 

Ar' = r 2 - r' = Ar + v 


L-l) Ar-y [(l-. a )Al^q 

a J v z a L J 


(56) 


For Ar jj v there is Ar'= a Ar (Lorentz contraction) while for ArJLy there is no contraction: 


Ar'* Ar 


2. Two events observed in the same point (Ar = r 2 - r x = 0) of the system S in the time interval 
A t *1 > appear in the system S' in the time interval A t' = /St/a (time dilatation), however 

now in the distance A r / = - v At/a = — v A t 


Thus, if Z 0 is the rest length of a rod, the length of the moved rod is 

/ = l Q \f\ — (v/c) 2 


(57) 


and, if r is the proper time of a clock which is at rest, the time of the moved clock is 


yjl — (v/ c) 2 

so that the moving clock will lag behind one at rest. 


(58) 
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D. Transformation of Particle Accelerations 


Differentiation of eqs. (13) and (14) gives the corresponding transformation equations for 
the acceleration: 


^ ^ ... • ?(-T^) 


dt' dt'/dt 


du dv/ dt' a 

a = = = a 2 


h) 


dt dt/dt' 


( l ^) 


Some special cases will be considered now. For a |[ u there follows 

a*v v 


a = a ' 


1 + a 


_ u»v \ 


and for a II u II v there is 


(59) 


(60) 


3) 


v 1 —UV/ C 2 


Setting a||v yields 


/ u*v\ a*v 

(‘-Tf) 


and setting aj_v (a»v= 0) gives 


✓ a 2 a 

n — 

_ ( « 

[1 — (u*v)/c 2 J 2 ~\1 -(u*v)/ 

Another case, ulv (u*v = 0) provides 


a '= a 2 a + 

a ‘ V (u V S 

c 2 \ 1 + a/ 


For the special case where v — ( v,o,o ) is in the direction of the positive x — axis the eqs. (59) 
and (60) yield 

= a 3 a x _ a 3 a' 

x (1 — u^v/c 2 ) 3 * (1 + u^v/c 2J 3 

a 2 tfy a 2 {a^v/cP)u y .. a 2 tfy a 2 (fl' i/c 2 j z/y 

a y =y ~ (1-u^v/c 2 ) 2 + (l-« x t/c2) 3 ; a y = y_ (l + a^/c^ 2 (l + «>/ C 2)3 


a 2 a z aH^a/c^a,, 
(l-^t/c 2 ) 2 ( i — «r t/c 2 ) 3 


a„=z = 


a 2 a 

(1+a; 


a/c 2 ,) 2 


a 2 («; a/c 2 ) a; 
(1 + a; a/c 2 ) 3 
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A very important case is obtained by putting u = v or u 0 (transformation to a rest system). 
This means that the primed system is centered in the moving vehicle itself. T he correspond ing 
time t' is called the proper time r . The quantities v = u and a ~ \]\ - (v/c) 2 - \J\ 1 - (u/c) 2 are 
now functions of the time t. The eqs. (59) and (60) are now 


a a + 


■ 2 l + o 


a a + 


(1 - a) u 


( 61 ) 


a = a 2 



a*'u u 

c 2 1 + a _ 


a 2 |a ^^-(1 - a) u 


( 62 ) 


and thus 


a*u 


(63) 


There are two special cases. Namely, for a || u or a'||u there is 

a ' = a/a 3 

and for a u or a / u there is 

a'= a/a 2 

In the section on dynamics there appears the function 


dt \jl-(u/c)V (1 -rf/c 2 )'/ 2 (l-uVc 2 // 2 C2 


(64) 


(65) 


(66) 


n a*u 9 a*u n 9 ^ 

a 2 a + — — u a 2 a + — — v 1 - a 2 ) u 

c 2 u 2 


where a = da /dt and u du/dt - a*da/dt = a*u has been used. Using eqs. (62) and (63) eq. 
(66) can be written 


d f u a 

dt \y/l —(u/ c) 2 ) a 


ou u / a*u u 

+ — = a a + 


.2 n 3 


c 2 1 + a 


(67) 


= aa' + {a'-a/a 2 ) = (l + a)a' — — «- 


For a || u or a '|| u there is 


, \ ) 

dt \\/l - ( u / c) 2 / 


a/a 3 = a ' 


while for aj_ u or a'J_u there is 


d ( u 

~dt ~ VVT - fa/^F 



= aa' 


( 68 ) 


(69) 
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2. RELATIVISTIC ROCKET DYNAMICS 


A, Definitions 

In accordance with the principles of relativity that the theorems of conservation of mass and 
momentum hold in all sets of coordinate systems in uniform relative motion (using the more 
complicated Lorentz transformations instead of the Galilean transformation equations) , it is 
necessary to modify the older Newtonian mechanics by assuming that the mass of the particle 
depends on its velocity. 

If m Q is the proper mass or rest mass of a particle moving with the velocity u the following 
definition equations hold. The mass is given by 


m = 


y 1 — u 2 - 7~c 3 


(70) 


This yields, for the momentum vector, 


p = mu 


m a U 

\Jl —u 2 Zc* 


(71) 


The total energy is introduced by 

E = me 2 = 


m n c * 


y I — u 2/c 2 yT — u^/c 2 


(72) 


where E q = m o c 2 is the rest energy. The difference between E and E o is the kinetic energy of 
the particle 


T = E ~ E o = = =7^ “ 2 (1+ 3/4 u 2 / c 2 +*••) 

The force is defined in relativistic mechanics as the momentum flow rate 

_ dp d . . d / 777 _ u \ 

F ~ dt = dt (mu) = ~dt \\/l - ui/cT~) (73) 

Tolman (Ref. 5, p. 46) says that ''the inclusion of m Q inside the bracket makes the expression 
applicable also in cases where the proper mass of the particle varies, as it might, for example, 
from an inflow of In 1934, he probably did not have in mind an application to fast-moving 

rockets expelling proper mass (exhaust gases). As in the author’s paper of 1955 (Ref. § ) where 
the Special Theory of Relativity was extended to systems with timely changeable rest masses 
(rockets), the proper or rest mass m Q will be treated as a function of the time t (but not of the 
flight velocity u) in this paper. Thus the last equation yields 


F - m 


du dm dm 

_ + __ u = wa+ _ ru = 


- . — 

VI - U2/C2 



± ( - 

dt \\/.l — u2/c2 


-) 


U 


(74) 


showing that in relativistic mechanics the force F and acceleration a will generally not be in 
the same direction, as it was in Newtonian mechanics. The power or work done on a particle per 
unit time will be defined here as the total energy flow rate 


dE 

c* — - 

d ( m 0 c 2 

V d ( 4 

dt ~ 

dt - 

dt \Vl-uVc2, 

J dt\sjr-ui/6i ) 


B. Transformations 


(75) 


The general Lorentz transformation (eqs. 1 and 2) provided the relations between the 
space and time coordinates of system S to those of system S'. The corresponding transformation 
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equations for velocity and acceleration were derived from them. It is very important to possess 
the transformation relations for certain other quantities of particle dynamics. From eq. (70) 
there follows 


m \jl-u 2 /c 2 - 772 ' yjl -u'Vc 2 = 77z o (invariant) 

Tlie propef time element (dr or dt Q ) is given by 

dt y/1 ~u % /c 2 - dt' yjl — u'Vc 2 = dt Q (invariant) 

A division of both equations shows the following simple relations 

dt dt ' 

m m 

By means of eq. (17) the eqs. (76) and (77) immediately give 

m dt jl — u ' 2 / c 2 a 1 + (u W)/c 2 _ / 1 -f (u • v)/ c 2 

Tfi dt' yl-u 2 /c 2 1 — fo*v)/c 2 a ^|l-(u*v)/c 2 

where 

a - \/l " y2 / c 2 

The transformation equations for these quantities of particle dynamics are, 

(1) for the mass: 

, l + (u'*v)/c 2 772 '+ (p'«v)/c 2 ,• l-(u*v)/c 2 772 — ( p* v) /c 2 

m — - m — = — — — • — - ; m = 772 — — — = — — 

a a a a 

(2) for the mass flow rates 

dm ^ dm' 1 + (u W)/c 2 m' /dv' v \ dm' „ ( o W)/ c 2 

dt ~ dt ' a + a c 2 / dt' + m 1 + (u '«v)/c 2 

or 

<&72 ' ^772 _ (a»v)/c 2 

IF ~ ~dt~ m T^TuTvIT^ 

thus 

, (° '•*)/<? _ (a»v)/ c 2 # a* v _ /m \ 3 

772 1 + (u *v)/c 2 ~ 772 1 - (u*v)/ d 2 ' a*v \772 / 

m ' u ' +m '^l (1 - a) iff + i] 

thus 


(3) for the momentum vector: 

772 a u'-f 772 V {(1 — -a)(ll W)/ Z 7 2 + 1 } 


P = 772 U = 


TTTu- v)7 c 2 


(invariant) 


or 


p =* .'+— r u j^ + n- p- + 4 

a L 2 c 2 J c 2 a 

a L v 2 c 2 J c 2 a 


( 76 ) 

(77) 

(78) 

(79) 


(80) 

(81) 

(82) 

(83) 


(84) 

(85) 
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(4) for the total energy: 


c- 2 , 9 1 + (u • v)/c 2 E + (p *v) nv E - (p*v) 

E = me 2 = m c 2 ; E - 

a a a 

( 5 ) for the force (momentum flow rate): 

vfFW dEn 

_ dp dp/dt' aF'+v{(l—a)(.F'-v)/v 2 +(dE'/dt')/c 2 } a ~ + c^U. +a + dt' J 
dt~ dt/dt' ~ l + (u'.v)/c 2 1 + (u'-vj/c 2 


( 86 ) 


(87) 


F' = 


a){ F*y)/v 2 — ( dE/dt)/ c 2 \ 

_ v TF»v dE~\ 

aF + ii\ 

1 - (u*v)/c 2 

1 - (u*v)/ c 2 


(6) for the power (total energy flow rate): 


dE _ dE/dt ' _ dE'/dt' FW) . (F.v) 

dt/dt' 1 + (u * y)/c~ 2 ' dt ' l-(u*v)/c 2 


( 88 ) 


(89) 


Instead of differentiating eq. ( 86 ) it is also possible to multiply eqs. (81) or (82) with c 2 
yielding 


dE dE' m'( a^v) 4 dE' _ dE m(a*v) ( 90 ) 

It "IF + 1 + (oW)/c2 ' IF ~ dt ~ l-(u.v)/c 2 

It is easy to see that the eqs. (89) and (90) are the same since 

F = -Jj fmu) = ra a + ^ u ; F ' = (mu') = m 'a' + u ' (91) 

The transformation eqs. (89) are not given in any standard book on Theory of Relativity, al- 
though they follow from the definition F = dp/dt or F = dp'/dt'. See, for example, the books 
of Tolman (Ref. 5) and Moller (Ref. 6 ) on Theory of Relativity. Instead of eq. ( 89 ) a further 
definition is introduced, namely 


dE/dt = (F*o) or dE'/dt '= (F'*ii') 

That would mean that in any system the change of kinetic or total energy per unit time is equal 
to the work done by the force per unit time. However, this is an overspecification. 

Identifying the system 5* ' with a rest system Sq (u ' = 0 ; u = v) the definitions dE/dt = (F*o) 
and dE'/dt - (F ' ) are correct only when the rest mass is not changeable with time. This 
assumption is always fulfilled when applying the Theory of Special Relativity to fast moving 
electrons, atoms or nuclear particles. However, this is not the case when applying it to fast 
moving rockets, expelling rest mass (exhaust gases). In the general case, these definitions are 
wrong and have to be replaced by the . correct transformation formula eq. (89) for the power. 

C. Transformation to a Rest System 

Putting u ' » 0 (u = v) in the before -mentioned transformation equations the following rela- 
tions will result when the prime is replaced by the subscript 0 (to designate that quantities 
belong to the rest system): 


m ^ E _ dt _ 1 1 (92) 

m E dt a a/1 - z/ 2 ?c 2 

000 v 
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and for the mass flow rate 


dm dm Q 
dt ~ dt n 


■+ m . 


fa 0 ‘ii\ dm Q dm m ^0*0^ 

V cV ; dt = Tt ^2 l TV 


( 93) 


thus 


The force is transformed according to 


a»u 


F + u 


P F u 

F ° = T~ + 


*u 1 dE „ 


(1 - a )-^- + 


2 dt 


,, v F*u 1 <s?E 

U - a) — = -t— 

u 2 c 2 a/ 


F U 

„ "l - « 


_F_ « 

a a 2 c‘ 


F *u zfE 

_Q _j_ £ 

1 +a dt Q j 

F > u _ dE_ 

1 + a dt 


( 94) 
( 95) 


thus 


(F*u) = (F. *u) + (1 - a 2 ) 


d E n I _ x 0 dm n 

-r -* = (F *u) + zz 2 -j-*- 

0 dt n 


( 96) 


The transformation law for the power is given by 


dE dE / x dE ( — x 

dT = dC +m ° = 


( 97) 


dE„ dE m , \ 1 

= t: o Ca-u) = o - 


dt 


dt 


dE 

dt 


- (F*u) 


( 98) 


The eqs. (96), (97), (98) were already given by the author in 1955 (Ref. 8). Writing eqs. (61) 
and (62) as 


a - a 2 


a a + 


1 + a 


a o »u 


= a A 


«. _d_ a )(^,u 


u f°' u \ 

l+a \c 2 J 


a a 4- (1 - a) ( ) u 


( 99) 


(100) 


and using these transformation laws for the acceleration or the proper acceleration* respectively, 
together with eqs. (92) and (93) in the definitive equations for the forces 


d , s dm - 

_^ u) = ma + — u ; F o = m o a o 


( 101 ) 


then the above-given transformation equations (94) and (95) are obtained. Dividing eqs. (97) 
and (98) by c 2 yields 


dm dm Q F 0 »u # dm n 
dt dt ' + c 2 * dt 


I f dm F*u 


a 2 dt 


( 102 ) 


The first of these equations can be obtained directly without using the transformation equations 
by differentiation: 


dm _ d frn\ m Q /a*u\ 1 

dt dt \a ) a 3 \cy + a 


1 dm. /a„*u\ dm„ 

dt ~ \c 2 ) + dt 


F *u\ dm. 
— ) 4 - — 2 - 


dt 
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Instead of splitting up the force, F, in one part proportional to the acceleration a and 
another part proportional to the velocity u (see eq. 101), F can be subdivided into one term 
due to the change of the rest mass with time and the rest term F* identical with F for parti- 
cles with constant rest mass. 


d (m n u\ d fu\ u dm 

~ di \^r) _77Z °5rw 


dm ^ 


h 14,11 L 

df- =F * + lt u 


(103) 


Due to eq. (67) 


d (U\ a /a*u\ u a /a *u, 

~dt W “ a + - a + \c*) u " aa o+ c 


*u / u 

2 \1 + CL 


thus 


and 


F* = 


t; Gi-^) " "• irtjr) " * (V9“ - 


F* = tzF +( = <zF„ + (l-a/^ 1 ^) u 


c V 1 + a 

Multiplying the last equation with u yields 

( c \ (c \ dE dE 

(F^u) ’ (F ° ,u) = Tt ~ l£ 

Putting eq. (106) into eq. (104) gives the equation of motion in the form. 

d\i c /F**u\ 

m -r-= F*- u 

dt * V c 2 / 


(105) 


(106) 


(107) 


Compared with the corresponding equation in the book of Miller (Ref. 6) eq. (107) has the 
additional term - a 2 (dm o /dt Q )v on the right side and the advantage that it is also applicable 
to systems with timely changeable rest masses (rockets). Introducing the force F, or the 
Newtonian force F 0 , into equation (107) yields [by means of eqs. (102), (105) and (106)]: 


and 






(1 — a) 



In the special case F q || u there is 


(108) 


(109) 


F* = F • F = F„ 


dm 
a 

dt o 


m a = a 2 F or a/a 3 = 


( 110 ) 


D. Application of Relativistic Dynamics to Rocket Propulsion 

In the following, the data for a rocket and its exhausted gases will be marked without a 
subscript in the system S of a stationary earth observer and with the subscript 0 in the rest 
system S Q of an astronaut centered in the moving rocket itself. Applying the transformation 
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eqs. (13) and (14) for particle velocities to the exhaust velocities v eo in the system 5' o 
where the thrust is generated and v e in the system S of an earth observer the following formu- 
las are obtained: 


v 


e 


av^f u {(1 ~ q)(u *^)/u 2 - 1} _ av^-f u{(u»v /c 2 V(l + a) — 1} 
l -(u'v^l/c 2 

qy e -m {(1 - g)(u*v e )/z^--f 1} _ av e + u Ku»v e (_ c 2 Ml + a) + 1} 

1 + (o*v e )/c 2 l+(u*v e )/c 2 


( 111 ) 

( 112 ) 


where the relative velocity between the two systems is u and a - \/l — u 2 /c 2 . Eq. (79) can 
be written, for this case, 


V I - t/Vc 2 " 1 a 1 + (u»v ft )/c 2 _ / 1 + (u*v e )/c 2 _ dm o (113) 

1 -^2 /c 2 ~~ 1 ~-(u*V eo )/c2 ~ a ^yj 1 -(u*V )/c 2 dm 


where dm and dm Q are the elements of mass flow after ejection in the two systems S and 
S Q , respectively. The element of mass flow before ejection 


dm yjl — v \/c 2 = dm Q \jl ~ /c 2 = dm* 
is invariant. The mass flow rate in the system S Q of the astronaut is 




dm n 

dt 


Vi - 


(114) 


(115) 


where /x* = dm*/dt Q is the mass flow rate in a system in rest relative to the exhaust gases 

of the rocket. The corresponding mass flow rate in the system S of the earth observer is 


P s - 


dm 


art 


1 - (u«v )/c 2 


* ^ w ’ft p / 

dt s/1 - v\ /c2 = /i ° VI - 


< 1 . 6 , 


in accordance to eq. (102) . The quantity 

Fo — f^o ^ eo 


VI -v*/c* 


(117) 


is the primary thrust force of the rocket in the rest system S Q where the thrust is generated. 
The corresponding thrust in system S which the earth observer would measure is 


F = /iv e = 


Vl-^ 2 /c 2 


= p; 


1 -• (u‘V M )/c 2 

V ; i - ^ Ve 


(118) 


or using eqs. (Ill) and (117) 


-«? 


av m + u }(1 - aKo'V^ )/a 2 - 1} 


Vi - V * A 2 


= aF 


+ u 


(1 -a) 




am 


in accordance with eq. (94). The results can be summarized in table 2. 


An important quantity is that part of the total energy of the exhaust gases which can be 
converted into useful work (kinetic energy). In the rest system S of the astronaut this conver- 
sion factor, first introduced by E. Saenger (Ref. 9), is given by 


(dE 0 /dt 0 )-(dE*/dt 0 ) 

dE 0 /dt 0 


- (**- = i - V 1 ~ v 2o /cl 


( 120 ) 
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and thus 


v w /c = y/1 - (1 - y) 2 = V/2 -y) 


Putting ot = 772 o /m into eq. (113) gives 


G?772 

772 



The equation of motion is 


d . d U ^772 - 

— ( 772U) = 772 -5 + U=F 

dt dt dt 


dm 


dt 


( 121 ) 


( 122 ) 


(123) 


or 


du / 

It = " (V e 


x dm 

+ u) IT 


thus, using eq. (122), 


du ~ - (v e + u) 


dm 

772 


= “ (v e + u) 



»'*eo'\ dm Q 
C 2 ) rn o 


Applying eq. (Ill) then yields 



Considering the special case || u equation (124) reduces to 


(124) 


dm _ 1 du 1 du 

m o v eo a2 v eo 1 ~'“ 2/c2 


(125) 


which can be easily integrated for constant exhaust velocity. Let be the rest mass of the 

rocket for t ~ 0(a = 0); the integration yields 


thus 



( 126 ) 


where r is the mass ratio of the rocket. The inverted formulas are 


u_ _ 1 - : (l/r) 2Veo/c _ r 2v ~ /c —1 
c ~ 1 + (1/r) 2 v ~ /c " r 2 v eo / c + 1 


(127) 


The eqs. (126) and (127) represent the fundamental relativistic rocket equatio n derived b y 
J. Acke ret (Re f. 10) in 1946 from the conservation law of momentum. Using a = y/l — u 2 /c 2 ; 
u/c = \J\ — a 2 eq. (126) can be transformed into 

l rn Q ( g 

r M q \1 + yjl — a 2 



(128) 
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This formula was found by R. Esnault-Pelterie (Ref. 11) in 1928 from a relativistic treatment 
of the special case of a rocket with constant proper acceleration a Q . He noted that this equation 
is valid for any law of mass consumption since it is independent of a 0 . In classical physics 
( c °o or u/c 0) the basic relativistic rocket equation becomes the well-known formula 


AT 

r = — — = e 

772 „ 


u/ v e 


or In r = 


since 


. . , . . u ln{l + u/ c)- In (1 — u/ c) 

limln r- lim — 7 

u/c-»0 U/C-+Q V ZU/C 


lim — 
u/c-»0 V eQ 


1/(1 + u/c)+ 1/(1 - u/c) 
2 



Representing 



1-:X 

1+ Xj 


c/2 v eo 


(04 y ^ 1 ) 


as function of x = u/c (0 ^ x 4 1) > a double differentiation shows that an inflection point 
(y " = 0) appears for 


The most favorable case = c (photons ) yie lds x. = V 2 and y. - 1 /yj 3 ; that is, the 

inflection point appears when u = c/2 and r = yj 3 

Before closing this part, it should be mentioned that in the special case F 0 || u 

F=F o - llo u=F 0 (l- ^ ( 129 ) 

H u - Veo 1S attained, then v e will vanish and also F = 0. From the technical point of view 
the thrust force F Q of the rocket in the rest system S 0 is naturally the matter of primary 
interest. 


xi = 


2v~ 


< 1 and 


/i = 


E. Motion and Mass Consumption of a Rocket 
Under a Constant Proper Acceleration 

The motion of a body under the action of a constant proper acceleration a Q in the direction 
of the velocity (a 0 || u), that is rectilinear uniformly accelerated motion, is known as hyperbolic 
motion. It was first considered by Minkowski (Ref. 12) in 1908, but was discussed more ful- 
ly by M. Born (Ref. 13) in 1909 and A. Sommerfeld (Ref. 14) in 1910. It was R. Esnault- 
Pelterie (Ref. 1) in 1928 who applied this kind of motion to rockets with constant thrust accele- 
ration and discovered the basic relativistic rocket equation for any law of mass consumption. 
Contributions to this problem were also made by Shepherd (Ref. 15), W. L. Bade (Ref. 16), 
E. Saenger (Ref. 17) and others. This presentation follows closely the author's treatment 
(Ref. 8). 

For a rocket in rectilinear motion with a constant proper acceleration in a system moving 
with the rocket (rest system) and thus a different one at each instant-eq. (67) can be written as 

It ( vi - U 2 / C r) = Ti-a 2 / C 2)3/2 = a ° (const - ) < 13 °) 
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or, integrated, 


y]\-U 2 /c 2 


-at (with u - 0 for t = 0) 


This gives for the velocity 


dx 

u = = c 


dt VI + ( a Q t/c) 2 


a n t 1 


1*3 A 


2-4 V 


5 t) 


and, thus, fora 

a= ^2- = V1-'« 2 /c2 = — = t, - T ] 

dt m v 0 £ yT+Ttf 

Eqs. (130) and (133) yield for the acceleration 


crgr-^KS**) 


rO')’ + &(H *-<»» 


“ " aoC(3 [1 + (a °/c) 2 j3/2 


« li- k»l(f^ a + 

° |_ 12 \c 7 1-2 2 2 \c 


-a t +••<• 


The distance is obtained by integration of eq. (132) with 


, a n 

1 + -2- X = 1+ -2- < 

c 2 c 


VI - a 2 /c 2 


Eqs. (135) or (136) give immediately 




-f *+l 

. C 2 


1 or (x + 


- (c f) 2 


Therefore, the world lines are hyperbolas in an x — t diagram and the motion is called hyper- 
bolic motion in comparison to parabolic motion (x - l / 2 a 0 t 2 ) in Newtonian mechanics. The 
inversion of the above-mentioned formulas yields for the time 


u/ c 

— i T 2 /c~ 


yi — a : 


i+ ^ x) = r^x ( 2 + ^ x 


The proper time of the astronaut follows from integration of eq. (133): 


t = — In - k / 1 + 

° C \J 


sinh r- 2* 


At) 

\c / 


“L 1 

The inverted formula reads 


_ j_ (* a .'s\ _k§_ A A 

2-3 \c 7 2-4.5 \c 7 


f°o , \ 

fc 

So to/c 



J 

e 



1 1 a 1 «„ 4 

= t n 1 + -z-r — 2 1 + -=-r — L + •« 

0 3 ! c ° 5/ c 0 


Inserting this relation into eqs. (133), (134), (132) and (135) yields the acceleration, velocity 
and distance of the rocket as a function of the proper or local time of the astronaut. 


2 4 

5 


cosh (a Q 
a o 

cosh 3 (a t /c) 


- 1 ' KrH T "* i: 

■■“"‘(f 1 ) 1_ 'T (?*'•) *TT G*') T 

(HH 


cosh 


a 2 

, 1 

fa Y 1 

/a \ 4 

= —ft t 
2 * 

1+ 3-4 1 

He V 3- 4- 5- 6 

H 


(141) 

(142) 

(143) 

(144) 


In addition to eq. (139) for the proper time t„ the inversion of the preceding formulas yields 

"(f) 




u / c 
u/ c 


- tanh 


= In 


1 + i? x + - x / 1+ (^ x ) T-1 =cosl,1 ( 1 + ^ x ) 


(145) 


A rocket travelling with constant proper acceleration a Q requires, in order to attain a given 
distance X, a time 




(146) 


in the system of classical physics; 

t ■■ 




2 + -S-X = TJ 1 + 


o„X 

2c 2 


= r 


l + 


1 + _M_ fSJL) 3 +... 

!/ 2-4 \ 2 c 2 / 2 * 4-6 \ 2 c 2 ) 


2 V2c 2 


(147) 


in the system of the stationary earth observer; 



O 


J cosh- Hi + fl„X/c2) 
V2 a o X77T 


= T 


1 


J_/£o*\ + _M_ f^ 2 LV_ 1 - 3-5 
2-3 \ 2 e 2 J 2 - 4-5 \ 2 c 2 / 2 - 4 - 6- 7 


sinh'Va n X/( 2 c 2 ) 

V« X/(2c 3T 

v O 



(148) 


in the system of the astronaut traveling with the rocket. In general t < T < t. Thus the astro- 
nauts undergo an immense gain in time from the relativistic principle. 

The preceding relations contained only kinematical data. The equation of motion in the 
system of the astronaut 


m a 


F =• 


dm Q 
dt _ 


or 




( 149 ) 


yields, after integration, the mass ratio of the rocket: 


AT 


... fs., ^77^] [i± y" — 


c/Foo 




c/veo 


(150) 


The thrust is given by 


F o = 772 = M 0 /r = AC <a e 

o o o o o o o 


a o*o/ v eo 


’ p o /m c 

The other kinematical data, expressed as functions of r, are: 


; a - \jl-u 2 /c 2 ~ 


2t 


f 2ve o/ c _ j 

f 2v 00 / c ^ l v f2veo/ c ^ 

.2 / , s ^2 


X- sl(L _ A _ iZL 

^ 0 \a 7 2<z 0 r 


v©o/c _ ^ j2 


v eo/ c 


‘° - ^ ln ' 


_ c / a/c\ _ c r _ 
' ~ \ a / ~ 2 0 r 


2 Veo/c 


- 1 


v eo / c 


(151) 


(152) 

(153) 

(154) 


All relations in this paragraph can be written in a dimensionless form by selecting c/a Q ' as the 
unit of time, c as the unit of velocity and cVa Q as the unit of length. With c « 3* 10 10 cm/sec 
and a Q ~ g o - 981 cm/sec 2 


c 2 
8 o 


9 - 1Q20 
981 


9.18-1017 cm 


9.18 - 10* 2 km = 6.14 - 70 4 A.U. = 0.97 light years 


c 


£o 


3 - 10 10 

981 


3.06- 10 7 sec -3(34.2 d ^ 0.97 years 

F. Motion and Mass Consumption of a Rocket 
Under a Constant Thrust (Constant Mass 
Flow Rate) 


In this case, first treated by the author (Ref. 8) in 1955, and also by Kooy (Ref. 18) the in- 
vestigation again proceeds from the equation of motion of the rocket in the system of the astro- 
naut traveling with the rocket, namely 


m n a _ 

(l _ u 2 /c 2 )3/2 = m o a o= Fo= V eo Vo (155) 

Taking a constant mass flow rate, fi o , the mass then decreases linearly with time t Q according 
to the law 

TO o = M„- Ho*o= «o (i- (156) 


Again taking the exhaust velocity v eo of the gases to be a constant, then the thrust, F 0 , also 
must be constant. On the other hand, F = F 0 (1 — w/v eo ) = (v^ — u) is variable. Combin- 

ing the two eqs. (155) and (156) the proper acceleration is 


(1 - u Vc 2 ) 3 ' 2 ~ a ° - m 


^° /M o 


1 - (no /M a ) t c 


(157) 
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or (since a - du/ dt = a du/ dt Q ) 


du 

1 — z/ 2 /c 2 


Hs/M o 

1 - (/x 0 /M 0 )t 0 




An integration yields 


c 1 + u/c . 

“o* T 7— = - V eo 

2 1 —u/c 60 


O-XH 


Thus, the reciprocal mass ratio is 


^0.^1^. iLo, ^ 
AL ° 


1 — u/c 
1 + u/c 


c/2 veo 


( 158 ) 


This is again the basic relativistic rocket equation which is valid for any law of mass consump- 
tion. The inversion of eq. (158) gives the velocity 


u = c 


1 — (1/ r) 2v e° /c 

1 + (l/r) 2v eo/c 


2 v eo / c q 

2 v e0 / c +1 


(159) 


therefore 


dt n rt 2(1/ r) Veo/c 

dt ~ V/i “ 1 +(l/r)2veo/o 


2 r Veo/c 
r 2 v eo / c q 


(160) 


In the two preceding equations, l/r can also be replaced by m 0 /M 0 or 1 - (/x 0 /M 0 ) t Q . The 
acceleration is easily found from a - a 0 a 3 using eq. (160). The time and the distance can be 
determined in the following manner, using eqs. (158), (159) and (160): 


J 0 a J ra 




1 -(l/r) 1 -( v eo/Q l-(l/ r ) 1+ ( Y eo/c) 
1 - (v / c) + 1 + (t> / c) 


for v ^ c 


_ J 

~ \ 


and 


x ~ c 


r 


M 

o 

2 /o 

X 


, 1 - (l/r) * 

In r + 

M n f 
= -j-2- Ini 

r 2 + (1 - l/r 2 ) 

2 

4 ^o 1 



for t/ = 


f o 

m/ c 


Me 


*„= c — 
cz / 1 


< 


c_ 

2 Ho 

s_ ik 

2 


1 -..(l/ r ) 2 - (v °° /c) 
_1- (feo /c) 

1 - (l/r) 2 


f u/c 

J ar* 


dr - 


1 -(l/r) 


1 + (Veo/c) 


1 + (^eo / c) 


for Veo ^ c 


In r — 


4 ^ 


In r 2 — (1 - l/r 2 )| for y = 


In these last equations, l/r can be substituted by m Q , or z/ according to eq. (158). 


( 161 ) 


( 162 ) 


(163) 

(164) 


In classical physics (c «? ; f /<»0) a = 1/ a - and 1 = For eqs. (159) 

and (163) the limit must be found by means of the Bernoulli-de THospital formula: 


— — = lim -L ^ i - -i- lim — = lim (2r 2€ In r) = In r 

Veo 6->0 £ r 2£ + 1 2 £ 0 t Z £->0 
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and 


x 


~£o * 


1-0/r) 


l-£ 


1-0 /r) 


i +e 


1-f 


l + i 


= lim 4- 
e-*o 2e 


2f - (1 + () (1/ r) 1 ' e +(l- f )(l /r) 1 + e 

l-€ 2 


= lim 

e-» o 


= lim 
(-> o 


1 - 


(1 + e)(l/ r) 1_f -(l- e )(l/r) 1+f 


2f 






•'-t' 


(f) «->eTK4-f 


= 1-J_+J_ln L = 1 _.JL fi_:in I s ) =1- L±i!Li 

r r r r \ r/ r 


Unfortunately it is not possible to express acceleration, velocity and distance also by the 
time t , as in the special case of constant proper acceleration, since eq. (161) cannot be 
inverted to give ^ as a function of t. Developing this equation in a power series 


t = 


2 ^o 


l - (l - ^< d /mJ‘- (V m/c) 

l - lv~ 7cT~ 


1 — (1 -a 0 t 0 /M) 1 + <’~/c> 


1 + U„/c) 



then the reversed power series reads 


t n = t 



11 — 9 




( 166 ) 


All preceding relations can be written in dimensionless form when choosing M 0 /^ 0 as the unit 
of time, c as the unit of velocity and c M 0 /^t o as the unit of length. The following three 
tables show numerically the characteristic difference between the case of constant thrust accel- 
eration a Q and the case of constant mass flow rate fi 0 (the exhaust velocity v co was taken 
constant in all cases). For the second case (fi Q = const.) the calculations were made for two 
special exhaust velocities, namely for v eo = c/10 and v^o = c (photon rocket). The assump- 
tion Veo /c = 0.1 is not typical for ion rockets; it is probably an upper limit. For fusion pro- 
cesses, the mass conversion is <0.009 corresponding to v eQ /c <0.134. 
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Data of rocket in free space 
without external forces 


In system S' of stationary earth 
observer 


In rest system S Q of astronaut 
centered in the moving rocket itself 



r k* 






Table 3. Constant Thjrust Acceleration a { 


P-c 

771 dt 

a= m = It 

a 

X 

V 


t 

t 

Vo/c 

T 

a o 

c%/a 0 

c /a Q 

c /a Q 

0.0 

1 

1 

0 

0 

0 

1 

1 

0.1 

0.9950 

0.9850 

0.0050 

0.1005 

0.1003 

0.9980 

1.1054 

0.2 

0.9798 

0.9406 

0.0206 

0.2041 

0.2030 

0.9944 

1.2247 

0.3 

0.9539 

0.8680 

0.0483 

0.3145 

0.3097 

0.9847 

1.3627 

0.4 

0.9165 

0.7699 

0.0911 

0.4364 

0.4240 

0.9715 

1.5274 

0.5 

0.8660 

0.6495 

0.1547 

0.5774 

0.5493 

0.9513 

1.7321 

0.6 

0.8000 

0.5120 

0.2500 

0.7500 

0.6932 

0.9242 

2.0000 

0.7 

0.7141 

0.3642 

0.4004 

0.9803 

0.8675 

0.8850 

2.3806 

0.8 

0.6000 

0.2160 

0.6667 

1.3333 

1.0986 

0.8242 

3.0000 

0.9 

0.4359 

0.0828 

1.2941 

2.0647 

1.4722 

0.7130 

4.3589 

0.995 

0.1 

0.001 

9.000 

9.950 

2.9932 

0.3008 

19.975 

0.99995 

0.01 

0.000001 

99.00 

100.00 

5.298 

0.053 

200.00 

1.0 

0 

0 

oo 

OO 

OO 

0 

OO 



Table 4. Constant Thrust v a /c =0.1 (Ion Rocket) 


= (x/c 

cM 0 / ^o 

0 

0 

! 

o 

T: X 

o 

t 

T 

o:o 

0 

0 

0 

1 

l 

0.1 

0.0269 

0.6334 

0.6329 

0.9992 

2.724 x 10° 

0.2 

0.0603 

0.8719 

0.8684 

0.9960 

7.599 x 10° 

0.3 

0.0816 

0.9604 

0.9547 

0.9941 

2.208 x 10 1 

0.4 

0.0930 

0.9935 

0.9855 

0.9919 

6.912 x 10 1 

0.5 

0.0981 

1.0051 

0.9959 

0.9908 

2.430 x 10 2 

0.6 

0.1002 

1.0088 

0.9990 

0.9903 

1.024 x 10 3 

0.7 

0.1008 

1.0098 

0.9999 

0.9902 

5.844 x 10 s 

0.8 

0.1010 

1.0101 

1.0000 

0.9900 

5.905 x 10 4 

0.9 

0.1010 

1.0101 

1.0000 

0.9900 

2.476 x 10 6 

0.995 

0.1010 

1.0101 

1.0000 

0.9900 

1.024 x 10 13 

1.0 

0.1010 

1.0101 

1.0000 

0.9900 

OO 


Table 5. Constant Thrust v Q /c = 1( Photon Rocket) 


0.9979 

0.9935 

0.9856 

0.9735 

0.9576 

0.9362 

0.9067 

0.8642 

0.7920 

0.5439 

0 


3 n/r 


t 

1 

J — Ll/ L. 

cM 0^0 

M 0^ 0 

O 

o 

0.0 

0 

0 

0 

0.1 

0.0047 

0.0 955 

0.0953 

0.2 

0.0180 

0.1847 

0.1835 

0.3 

0.0394 

0.2701 

0.2662 

0.4 

0.0690 

0.3547 

0.3453 

0.5 

0.1080 

0.4414 

0.4227 

0.6 

0.1591 

0.5341 

0.5000 

0.7 

0.2278 

0.6396 

0.5799 

0.8 

0.3271 

0.7715 

0.6667 

0.9 

0.4993 

0.9730 

0.7706 

0.995 

1.2476- 

1.7464 

0.9499 

1.0 

OO 

OO 

1 


1.1054 

1.2247 

1.3627 

1.5274 

1.7321 

2.0000 

2.3806 

3.0000 

4.3589 

19.9750 


PART II 


RELATIVISTIC PERTURBATION THEORY OF AN ARTIFICIAL SATELLITE 
IN AN ARBITRARY ORBIT ABOUT THE ROTATING OBLATEO EARTH SPHEROID 
AND THE TIME DILATATION EFFECT FOR THIS SATELLITE 


SUMMARY 


In this part Einstein’s general theory of relativity (gravitational theory) is applied to the 
motion of an artificial satellite revolving in an arbitrary orbit around a central body and the time 
dilatation effect for this satellite is given. This relativistic perturbation theory is based on 
Einstein’s general field theory, differential geometry of non-Euclidean spaces, potential theory, 
and celestial mechanics. The short periodic perturbations are excluded by using time average 
values over a revolution. The secular and long-periodic (non-relativistic as well as relativistic) 
perturbations of the osculating orbital elements, which represent deviations from the elliptic 
orbit, are presented here for the case of a rotating, non-homogeneous, oblated spheroidal central 
body. This is an extension of the work of Einstein (1915) who considered motion around a mass 
point as well as the work of deSitter (1916) and, independently, of Lense and Thirring (1918), who 
treated the relativistic motion around a rotating, homogeneous, spherical central body, omitting 
the terms due to the square of the angular velocity. 

A formula for the relative difference of the time rates of a satellite clock, compared against 
a standard earth clock (time dilatation effect) is derived for orbits of any eccentricity and equa- 
torial inclinations, thus extending the paper of Winterberg (1955), Singer (195b) and Hoffmann 
(1957). 
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3. APPLICATION OF THE GENERAL THEORY OF RELATIVITY 
TO ARTIFICIAL SATELLITES 


A. Relativistic Perturbation Theory 

In this section Einstein’s general theory of relativity will be applied to determine the 
motion of an artificial satellite revolving around the rotating earth as well as the difference in 
time rates of a satellite clock and a standard earth clock. 

In Einstein’s general theory of relativity gravitation is determined by the 10 different com- 
ponents of a symmetric covariant tensor of the second rank called the fundamental or metric 

tensor. These components g kJ are functions of the coordinates x t ; x 2 ; x 3 , x 4 and they appear 
in the formula for the four dimensional line-element of the non-Enclidean time-space world, namely 

ds2 = | r 8k ‘ dXk dx ' (gk > = 8 ‘ k) (167) 

In the following, a spherical polar coordinate system (x t = r; x 2 = 0; x 3 = <fi; x 4 ~ct) will 
be used, where r is the radius vector, 6 the longitude, c£> the latitude, t the time and c the 
light velocity. The line element has the following form 


ds 2 - g u dx 2 + g 22 dx* + g 33 dx 2 3 + g 44 dx 2 4 + 2 g. 24 dx 2 dx 4 

=- (1 + a) dr 2 - ( 1 + jd) (r 2 cos 2 0 dQ 2 + r 2 d tfi 2 ) + (1 + y)c 2 dt 2 + 2g 24 d0 cdt 
- “ (1 + fi)do 2 + (fi - a) dr 2 + (1 + y) c 2 dt 2 +2 dQ cdt 


thus 


g u = - (1 + a) ; g 22 = -(1 + 13) r 2 cos 2 cf} ; g 33 = -(! + / 3)r 2 ; g u = 1 +y 


_ , . . , b cos 2 

g ik = 0(1 4 k) except g 24 = g 42 = * 


(b - const) 


(169) 


where a and are functions of r alone while g will be assumed to be a function of r 

and 0 because it is needed to the second order. 

The fundamental metric tensor has now the following covariant components: 


/ &11 &12 &13 &14 \ 

1 


1 8>n 

0 0 

0 \ 

&21 8 22 8 23 &24 


0 

822 0 

824 ' 

& 3 t 8 32 &33 834 

, J 


0 

® 833 

0 

\ 841 842 843 844 / 

1 

\° 

842 0 

844 


30 



with the determinant 


8 =\g kl \ = gi 


> 11 


£22 ^ 824 


833 0 


0 824 

® 833 ® 

8 11 8 22 


+ 8 n 842 


842 ® 844 


0 g 44 


833 0 


8 11 &22 &33 8 44 &11 &33 &24 ~ &11 &22 &33 & 44 

The contravariant components are now the minors of the correspondent covariant components 
divided by the determinant, namely 

ki = minor of g kJ cofactor of g kt 
8 8 


Therefore 


11 1 

8 = 7 


822 0 g 2 4 

0 8 33 0 

&42 ® 844 


1 33 1 

T u ' 8 -T 


811 0 0 

0 822 824 

0 8 42 844 


_ 8j 1(822844 824 ) __ ^ 

g U g 33 ( ^2 g 44- g 24 } *i 


y 33 




*« 0 0 
0 s 33 0 

0 0 g 


44 


- MU 833 8 44 


811 8 22 8 33 844 8 tl 8 33 8 24 


%22 824 / 844 


& 22 


o 44 = — 
8 8 


811 0 0 

0 *22 0 

0 0 g 


33 


81 1 822 833 


1 


811 822 833 844 8 tl 8 33 8 2 4 844 8 2 4 ' 8 2 2 


g 24 m g 42 = L 


8ll 0 0 

0 0 824 

0 S 33 0 


Mu 8 .3 3 82 4 824 _ _ 824 

811 8 33 ( 822 844 “* 824 ^ 8 2 2 844 “■ 824 82284 


The Christoffel three-index symbols of first and second kind can now be defined as: 

r _ \ kl } = L(§Ju^ + ifiL.i^VPl 

,kt L i J 2 \dx t + dx k dx t ) "L* J 


( 170 ) 


and 


r f 

ki 


■ f/} - Jy M- * Hr - {?} 
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32 


r*3 l d £>n , _ -p3 _ 1 d g 

^ i ' r ■*■ ^ i 


r 3 = 

22 


2 g 33 d i 


_ J _ ±822. ■ 


11 


13 31 2 g 33 d x t 




; r J = r J «- 

2 -# 42 


i <9g 




2 £33^*3 


p 3 _ A ^^33 .p^ = _ d g 

35 2 g 33 <?* 5 ' 44 2 g J3 d * 3 

and the following Christoffel symbols vanish 


p 4 _ _ L (± &24 _ 824 d 82 2 \ ; p * = p < = J , dg 44 

12 21 2g 44 \d x t g 22 dx 1 J 9 14 41 2g 44 dx j 

r< = r< = -?_ §24 dgil) 

23 32 2 g.. \dx, g „ <9*,/ 


p 4 _ p * _ _i d g dA 


34 


43 2 g 44 d l 


r 1 = 

r‘ = 

r' = 

r' = 

II 

T_ 

r* = 

r‘ 

= r J 

= 0 

l 2 

1 4 

21 

2 3 

32 

34 

41 

43 


"1 

II 

r 1 = 

r' = 

r‘ = 

r 1 = 

r 1 = 

r' 

II 

= O' 

1 1 

1 3 

22 

24 

31 

3 3 

42 

44 



By contraction of the Riemann-Christoffel curvature tensor of rank 4 for the non-Euclidean 
space-time world a symmetrical curvature tensor of rank 2 is obtained by the explicit expres- 


krl - R kl - R lk 


!>■ 

r 

Using the relation 


y r - r - ..y^y y r - r . 

^ rs Z j d X f / * * r 


(172) 


2 r : ■ 421 




r r r 

the contracted curvature tensor can now be expressed by 

iiwur 


d In |g„ | _ 1_ (9 In |g| _ d In V 

k ~ d x k 


- = iy 

2£j dx u 2 dx u dx 


kt Ik d x d 


TjT_ V r ' a W|gl _ Vr s r r 

X, 2 j ki dx t 2 jdx r Zj kr Is 


(173) 


_ <? 2 lnVlg| f 1 * <?lnVlgl t.3 ainVIfl 31^, <9r, 

— i .. n „ 1 n 1 ^ ~ — — -z — 1 


d x k d x t 


kl 


d x t 


kl d x. 


d x t 


d x. 


hL. 


+ [r' 

r 1 + r 2 r* + 

r 3 r* + r 4 r 1 1 

+ r r 1 

r 2 

+ F 2 

r 2 

+ r 3 

r 2 

+ 


r 

L kl 

11 kl 12 

kl 13 kl 14 J 

L k2 

u 

k2 

12 

k2 

13 


k2 


+ [r 1 

r 3 + r 2 r 3 

+ r 3 r J +r'r 3 1 

+ It 1 

r 4 

+ r 2 

r 4 

+ r 3 

r 4 

+ 

r 4 

r 

L k3 

It k3 1 2 

k3 13 k3 14 J 

l ** 

1 1 

k4 

1 2 

k4 

1 3 


k4 
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and thus for l = k 


<3 2 WH r j ^ ln Vlsl r 3 d inV|g| _ d rl k d r/t 


k/c d x 2 

k 


kk d X . 


kk d X , 


<9 x^ 9x ; 


2 

k3 


+ (rM 2 + 2 r 2 r 1 +2 r 3 r 1 + 2 r 4 r 3 + (r 2 ) 2 +2 r 3 r 

V kl) kl k2 kl k3 kl k4 \ k2 J k2 1 

+ 2 r 4 r 2 + fr 3 J 2 +2 r 4 r 3 + (r 4 ) 2 

k2 k 4 \ kj k3 k4 \ k4/ 

1 3 

Neglecting higher order terms the ten components of the contracted curvature tensor are (T = 


= r 3 = r 4 = o): 

33 23 


it dx 


33 d x„ 


pi d In \/|g| 

+ (r J ) 2 

+ (r 2 f +(r 3 ) 2 

li d x 

dx n 

12 13 

p3 <?W|g| 

d r 22 d r 22 

+ 2 r 2 r 1 + 2 r 3 

22 d X 3 

d x t d x 3 

21 22 22 

pi <?W|g| 

2 r 3 r 1 

+ (r 2 ) 2 + (r 4 ) ; 

33 d X 

dx t 31 33 

32 34 


14 


R = _ r 1 <?Wlgl _ r ^ i_WJgj_ _ d i \ l 4 _ d r 4 3 + 2r * r t + 2 r 4 r 

44 a a r) r A A d X - 9 X. f 9 X. 41 4 4 43 


44 <9x ; * 44 dx 3 dx t 

R i3= R 3i= - r 3 --- ln /l g l + r 2 r 2 + r 2 r 4 + r 4 r 4 

OX t dx 3 13 dx 3 12 32 14 32 14 34 


3 

43 44 


R 24T R 42 


_ r * dWlsl r 3 dinVM _ $ r 2 4 _ , r 2 r 1 +r 4 r 1 


dx t 


24 dx 


$x t d x 3 


21 42 21 44 


+ r*r 2 +r 5 r 2 + r 2 r 3 +r* r 4 +r 3 r 4 

22 41 22 43 23 42 24 41 24 43 


R 12^ R 21 0 * R 1 4~~ R 41 ~~ 0 » R 23~ R 32~ 0 ' R 34~ R 43~ ® 

The most important problem in Einstein’s theory was to set up the general equations deter- 
mining the gravitational field variables or the g kl when the distribution of mass is given- 
Einstein finally solved this problem in 1915 after several attempts and succeeded in finding the 
general field equations in covariant form corresponding to Poisson’s equationin Newton’s 
mechanics. Based on the theorems of conservation of energy and momentum fie (Ref.l 9) found 

R ki~~ 2 R & ki = ~ K T kJ (*174) 
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where 


( 175 ) 


R ■ EE 8 *' * 

k 1 

is the curvature scalar and 

Tkt = gtk8st ^s~ (176) 

T S 

is Einstein’s energy-momentum tensor neglecting the small contributions of pressure and elastic 
stresses ( p = density). The constant of proportionality 


k - 


8 nG 
c 2 


(177) 


follows from a comparison of the general field equations with Poisson’s equation ( G is Newton’^ 
gravitational constant). 

Multiplying Einstein’s field equations (174) by g kl and summing over k and / (con- 
traction) then follows 


R ~ k T 


because 

r = YiYi gkI Tkt and X 8kI = 4 

k 1 k 1 

Therefore, the field equations can be also written in the form 


r m = ~ k ( T ki 


1 

2 &kl 


T) 


(178) 


(179) 


( 180 ) 


Neglecting the pressure gives 

T U= T 22- T 33= 0 ■ T 44 = p (1 +y) , T = p (181) 

Einstein’s field equations are a system of nonlinear partial differential equations of the second 
order which must be solved simultaneously to obtain the components g kl of the metric tensor. 
An approximate solution for a field with spherical symmetry was first given by A. Einstein 
(Ref. 20) in 1915, making the assumption that 


V |g| =1 or /3 + - (a + y)n=0; 

2 m 

a = - y , /3 = 0 , y = - 


where 


G M = f 1.475 km (Sun) 
c2 [ 0.4435 cm (Earth) 


(182) 


is the gravitational radius of th e central mass M. In 1916, K. Schwarzschild (Ref. 21) gave 
the correct solution (with \f | g| - 1). 
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i In* 2m 

1 + a = ; = 0, y 

1 + y ^ 7 r 

The line element has now the form 

ds 2 =- _^ 7 - — r 2 (cos 2 cfj d0 2 +d<j> 2 ) + fl- c 2 dt 2 ( 183 ) 
Introducing harmonic or isotropic coordinates r , 0, <£, c£ defined by the transformation 

r =7* ^1 + « r" + m (^84) 


the line element assumes the form 

(l + ^j 4 (dF 2 + T 2 COS 20 d 02 + 72 dc/> 2 ) + 0 ' C* dt* (185) 

Including only first order terms in a - /3 and second-order terms in y and neglecting the small 
difference between r and T the line element is now given by 

ds 2 ~ ~ (l + (dr 2 + r 2 cos 2 <f> d 6 2 + r 2 d<fi 2 )+ ^1 - 2 ^ c 2 d t 2 (186) 


thus 


a = j6 = -y 


y = ~ 


2 m 

+ 


T 


2m 2 
7-2 


This particular form was derived and used by deSitter (Ref. 22) in 1916. It determines the 
metric of the non-Euclidean space-time world in the neighborhood of the mass M and the gravi- 
tational field and thus governs also the motion of satellites around this mass. 


In 1916 W. deSitter (Ref. 22) has also shown how the motion of satellites is influenced by 
the rotation of the central body according to Einstein’s gravitational theory by also introducing 
the component 


T 24 ~ ~ P ^ r 2 cos 2 cf> 

of the energy-momentum tensor in his field equations (12 = angular velocity of the central body). 
The solution then gives the component 

b cos 2 <£ 

g 24 = — — (b= const.) 

In 1918, J. Kense and H. Thirring (Ref. 23) independently solved the same problem using 
Einstein’s linear approximation solution for weak fields (Ref. 24) and they arrived at the same 
result as de Sitter. All these authors neglect terms proportional to 12 2 . 

In this paper the g kJ will be determined including terms due to 12 and 12 2 as well as 

terms due to the oblateness of the central body. This then leads to the relativistic perturbation 
theory for the motion of a satellite around a rotating oblated central body. Again, Einstein’s 
linear approximation method for weak fields will be used. 
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A Cartesian coordinate system (xj = x; x 2 = y; x 3 = z) with x 4 - ict for the imaginary 
time coordinate will be used. In this system all the g^ have the same value — 1 in zero 
approximation. If 

&ki ~ ^ ki + y ki * ® ki “ Q ( jl i ) (187) 

and 

1 4 

V kl = y k , + 2 ^ y aa ( 188 ) 

a =i 

then Einstein’s solution for weak fields is given by 


kl 



Tjg dV„ 
A 


along with the energy -momentum tensor 



ill. 


l kl 




= 9 



(189) 


( 190 ) 


the volume element 


d V 




dx' dy' dz' = i 



r' 2 dr' cos <f> ' d <f>' d 0' 


(191) 


and the distance of the mass element dm' from the attracted point (x - r cos <£ cos 0; y= r cos^# 
sin 0 ; z = r sin <f>) under consideration 

A = [{x - x') 2 +(y - y') 2 + (z - z') 2 ] 1/2 = r [ 1 - 2 cos a 4- 2 ] 1/2 (192) 


with 


cos a ~ sin ' sin + cos ' cos <£ cos (O' — 6) 
From potential theory for r> r' 


i = f 2(f) p * (cosa) 


(193) 


(194) 


where, according to Legendre, the surface harmonic P n (cos a) is given by the addition theorem 

n 

P n (cos a) = P„ (Sin 0 ) P n (sin $') + 2 ^” ~ ^ 7 K < sin< £> P n ( si « <f> ' ) cos s (6 d) 

(195) 


a=l 


and 


s , £3? a P (sin<£) 

P (sin <p) — cos® <£ — ■■ n— - , -L- 


£/(sin<£) * 


( 196 ) 


are the associated Legendre functions of Ferrer* 
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Under the assumption that the central body rotates around the z- axis the velocity components 
of the mass element are: 


dx } dx . Oy .Or 

= — z — =— = z — = 1 cos 0 sin 6/ 

0x 4 cat c c r 

dx . dy' . Ox' .Or' 

— 5 — * = — 1 — r = — l = — 1 COS 0 cos d 

dx 4 c dt c c ^ 


dx j q t ^ ^ 

dx 4 


dx. 
4 

dx 


The components of the energy-momentum tensor are therefore 

- (~~) cos 2 <£' sin 2 O' cos 2 $ sin 0'c os 0' 

cos 2 4>' s ^ n 0 * cos 0 ' - cos 2 $ ' cos 2 ® ' 


Tki ~ P o' 


c 


K ds J 


z^^-lcos <56 'sin 6 ' 


./0 r' 


-J cos ' cos 0 ' 


Introducing eqs. (190), (191) and (194) into eq. (189) yields 


y ki 


4UV C • / dx A 5 dx k dx /rV , 

—r /, | ‘{if) lf 4 iT 4 (t) P n( C °“’> d « 

jrlmafinn flip re is 


11 — u u 

In first approximation there is 


ds 

dx 4 


dxj 

ds 


dx 4 3 
ds 


thus 


CXJ M 

, 4G\^ I dx k dx 1 /r"\ n n , , , 

y« - 77 / , ) 3 ^ 3 ^( 7 ) 

n = 0 ^ 


0 i 


(^-) cos 


sin (T 


0 - I^jcos 0 ' cos 0' 


0 

0 


1 

(197) 


(198) 


where the components of the tensor (dx k /dx 4 )(dx j /dx 4 ) are given by the matrix ineq. (197). 
It follows at once that 

Yi' 3 = Y31 = 0 ’ Y23 = Y32 = 0 Y33 = 0 •' Yu = Y43 

while the other components are given by 


' = 0 


Yn =- ^ (f) 2 ^ 7»£ r/n + 2c ° S ^ / P " (C ° Sa) ^' 

n~ 0 


(199) 
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<N 

0 

II 

/ 

4 G 


c 2 r 



714 

24 

= ± i 


® St-/ "” + 

midp- 

n — n *•' n 


(cos a) ^772 " 

2 cos 2 cf) ' sin Q ' cos <9 ' P n (cos <7) ^772 " 


/ sin Q n , sis 
cos 0 AS P„ (cos cr) am 

^ cos# n 


(200) 

( 201 ) 

(202) 


n ~ 0 ~ 0 

The expression for P n (cos cr) in eq. (195) and the mass element, dm' ~ p r' 2 cos <fi' dr' d<f>' d0', 
will now be introduced into eqs. (197) to (202). The integration with respect to the local radius r, 
will be made from 0 to R ; with respect to the latitude, 0 ", from - n/2 to n/2; and with respect to 
the longitude, Q \ from 0 to 2 tt. For a spherically symmetric model the density is independent 
of the longitude, Q\ and will be assumed to be given by p = p (r'<£'). Therefore the integration 
with respect to Q' can be performed immediately. Ineq. (199) there occurs the integrals: 


/ 


2 77 

sin Q' 
cos 0 ' 


d0'= 0 


/ 


2 77 - 

sin 2 Q ' 
cos 2 0 


, cos S ( 0' — d ) dQ 


+ 4 

77 

+ T 


i/. 


■if: 


[1 qp cos 2(9'] cos s (#" — Q) dQ' 


{ cos [ (s + 2) 0 s #] -f cos [ (s - 2)0' — s(9j| dQ' = 
cos 2 Q (s = 2) 

0 (st 0 ) 


Therefore 


T11 

22 


4 G 


00 

a 


4 

(“)12 

l n = 0 */ 0 

/ 


(n - 2) / Pg (sin <ft) 

(» + 2) / 2 rn 


cos 2# 


+ 2 cos 2 qS" P ^ (sin <f>')dm 
cos <fi' pi (sin ^ 


(203) 


r ' n + 2 


0 dm jj 


The following integral appears in eq. (200 ) : 


thus 


cos s (Q'~ Q) dQ'= 0 


f 

0 

i- st f 

n = 0 ^0 


^44 = 

In eq. (201) the following integrals can be found; 

> 2 77 - 

sin Q ' cos Q' dQ' - 0 


dm 1 


(204) 


* 


27T 


r 

# sir 

4/0 


sin (9 "cos (9 "cos s ^ % 


u 

1 /T- 77 
= y 7 sin 


/ 


sin 2(9" cos s (#"—#) ^(9" 


s m [(5+2)9 - s 6 1 ] - sin [(s - 2) 0' - s Q] J d Q 


77 

, f 2 sir 

l 0 


2 sin 2 0 ,(s = 2) 
(^2) 
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4 G /f 2\ 2 


therefore 


In eq. ( 202 ) occurs the integrals 



0 ? — 2 ). P n (sin 0 ) s ] n 2 r' 1 + 2 cos 2 P £ (sin <fi')dm' 

\n + 2)/ 2 r n I 


f27T 

sin ( 9 ' 


, ^<9' = 0 


sin ( 9 / 


, COS 5 (0'- 0) <*0' = 


■ «jf (- f 


sin 0 . 

n* fs=l) 


0 (s =# 1) 


(s + 1 ) e'-so - 


|7s - 1) 0'- sfl] J 


, . 4 g /o 

y l 4 “ ± 1 Vc 


O O 

E 


(72 — 1 ) / P* (sin <£) sin (9 
(72 + 2 ) / r cos 0 


r n -r 1 cos 


cf)' P 1 (sin <f>')dm' 


Introducing 


n M P 


[ R n + 2 I 


r'" + 2 cos 2 <^' P n (sin <j}')dm' = 




4 - -35 '”^ 2 ' 


i_r 

W P n I 


r ' n P n (sin 0 ') afo / 


-i f«=o; 

0 (72=1) 


n MR«+1 


' n + 1 cos (ft) ' P 1 (sin / ) ^772 / 


0 

_l r.. 

R n + 2 I 


2 cos 2 <£' P 2 (sin<£') dm' = { 


ST» - r -| '”’ 1; 

= 0 (72 = 2) 


(72= 2) 


0 

(approximate values refer to a homogeneous body) and using the gravitational radius m = GM/c eqs. 
(203) to (206) become 


= 


4 772 A]P\ 2 


/ /P\ n 


KY‘ . . ,* ( 72 - 2 )/ L /R\ n i 0 

r) p n (sin *> + (^T 277 ^ (t) p » (sin ^ cos 20 


, 4m /flRN 2 ^ 
«i2 = yi2= — (— ) 


^ K (yj P n 


Srlj? t (fJ p - (sin # sin 20 
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*14“ *14 = ± Z — — 


14 ' 14 

24 24 


and therefore 

4 


^2 


sin 0 
cos ( 9 .J 


( 210 ) 


n= 1 


2 ^ ['-2 ©" p " <>i "" , - c^2 '^ r p " <sin ' ) 

O = 1 n= 2 n= 0 


a= 1 

Because 


there is 


8kk ^ + X/cJc 9 


*2 


2 Z Xaa 
a = l 


g \\ ~ ~ 1 “ ^ 1 “ 2 Jn 0 ) P " ^ ^ + (x) (n + 2 )/ L n( 7 ) P n (sin ^ ) C ° S 26 


n= 2. 


«33=- 1 - ~ T 1 - p (sin ^° 2 /n (?) p " (sin< ^ ,) 

L u=2 n-'Q 

s '.(*r 


(sin < 56 ) 

\* ' \ ^ / / 1 \' / 

^ = 2 n = 0 

Taking only the terms up to n =2 then the components of the fundamental tensor are 


g ll - - 1 - “7^ |^ 1 -/ 2 (“) P 2 ( sin <^) + (~t) ^ (~) COS 2 <£ cos 2 0 

g 33 = - 1 - ^ { 1 - J 2 ( 7 / P 2 (sin <£) - (^) 2 r + A 0-J P 2 (sin<£)J J 

( 7 / p 2 < sin <£>)] 


- 1 + 


o , • JS I ft K N 2 
P 2 (sin <£ ) + ( — — 


«12 (^“) (f) cos 2< f> s i n 2 d 

2m /QR\ /rs si 

I 1 — I cos (h 

r \ c / \r/ ^ cc 


r+A 


+ . 2 m ^QR 
_ 2 


sin 6 
cos 0 


(211) 

( 212 ) 

(213) 

(214) 

(215) 

( 216 ) 

(217) 

(218) 


In order to have a spherically symmetric field the very small terms proportional to m/r and (QR/c) 2 
will be neglected. In g 44 only, the term will be retained and the second order term 


2m 2 / r 2 - [1-2 J 2 (R/r) 2 ? 2 (sin<£)] 
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will be added according to de Sitter (eq. 186). The reason why g 44 is required with higher accuracy, 
is that it appears in the equations of motion multiplied with the large factor c 2 . The line element can 
now be written 


ds 2 = g u ( dx 2 + dy 2 + dz 2 ) - g 44 c 2 dt 2 + g 14 dx icdt + g 24 dy icdt 

= - Jl+ - (7 / ^2(7)]} (dx^dy^dz 2 ) 


+ 2 (”) l ~ 2 h(jr) P 2 ( sin <f>)yc 2 dt 2 

- — (— ) r (—) d x cdt+ — - d y cdt 

r \ c / \ r / r r \ c J \r / r 


Because 


da 2 = dx 2 + dy 2 + dz 2 = dr 2 + r 2 (cos 2 ft dO 2 + dft 2 ) 
x dy — y dx = r 2 cos 2 ft dO , — sin ft 

the line element can be written 
ds 2 = 


= -{1 + ~ [1-72(7/ p * (sin ^]] {dr 2 + t 2 cos 2 ft dO 2 + t 2 d ft 2 ) 
/l- 7 ^- j 1 - /2 (7) P 2 ( S 2 a &+(ftft) [r + A ( 7 J P 2 ( sin < 7 >)]] 

2 (/9 [ i “ 2/2 0 Pa ( sin ^)] V 2 * 2 


^ 2m T fttR 

T 


R cos 2 <f> dO c dt 


( 219 ) 


( 220 ) 

( 221 ) 


( 222 ) 


This is now the square of the four-dimensional line-element for a non-Euclidean space outside a rota- 
ting oblated urhomogenous spheroid of mass M. Writing 

ds2 =&i dx l + S 22 dx 2 2 +g 33 dx 2 3 + g 44 dx 2 +g 24 dx 2 dx 4 

where now 


x 1 = r ; x 2 = 0 ; x ft ; x 4 = c t 


then the components of the fundamental tensor are 


Su =“ j 1 + ~ 1 - / 2 (ft) p 2 Csin^)j] ; g 22 = gjj r 2 cos 2 <f> 

2 m y, / £1 R\ 2 

S 3 3=«xi r ’ S 24 = g 42 = — r -j R CO S 2 ^ 


5 44 ^ 


P 2 (sin ft) + 




P 2 ( s i° 


+ 2 


(ft) 1 - 2 7 2 (“) 2 p 2 <•**) 
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Their derivatives are given by 


and 


dg n _ 2 m 
dr r 2 


1 - 3 J 2 (j) 2 P 2 (sin 0) 


dg. 


T dgj 


= 2 T cos2< £ Un + 2 57 11 ) ; 


1 d_ 
§ 22 d 


§22 = £ 


Su dr 


d &24 
a r 


r f— ) R COS 2 ^ = _ ^24 


v c y 


d g 33 

(9 r 


= 2r(g 11+ ' fill.) ; J_ |4»_ 2 J_|ln 

2 / <5 33 r 5u 


dg 

(9 


= 2w f 

r - r’[ 


1-3/ 




2 v r y p 2 (sin + 



+ 3A (-) P 2 (sin </>) 


- I 1 ' 4 ' 2 


P 2 (sin <£) 


ill 


^ g 1 1 ^ ^ r f n1/_* / \ u T / f\ 1 

= 7~ ^Vjy p 2 (sjn $ = T~ /a (j) sm <?> cos <f> 


6 m T /R 


§_ g ; 


<9 0 


22= 2 r 2 sin </> cos </> f- g u + -1- cot <j> |^ u ') ; — — j^ 22 - 


d<£ J ’ «22 d</> 



__ 

1 

dgu. 

— 2 tan 





£11 

<9 </> 



d g 2 4 

(9 

I 

II 

r ( 

'Q, R\ 
c ) 

R sin $ cos <56 

1 = - 2 g 24 tan $ = 

2 r tan <f) 


= r 2 ll 

Ll ; 

1 

. ^ = J_ 

. djLu. 


d<f> 

d <£ 


£ 3 ; 

i d </> g n <?<£ 


dg 44 
<?<£ 

11 

1 ^ 
I s 

h 

_ 

C 1 - 

^ c m;J 

| /P\ 2 • , 

1 — ■ J sin <p cos 0 


It is now easy to write the values for the Christoffel symbols of second kind: 


r 1 = — 

u 2 g 


li 


dgn ^ 1 dg n m 

dr 2 dr r 2 


1 ~ 3 h (f) p 2 ( sin <t > ) 
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p 1 _ _i - ill- = _ — J (—') sin d> cos 

13 2 g u d<f) 2 d r 2 \rj 


p 1 = _d ^Al2 

22 2 g u dr 


T COS J 


<k ( 1 + ^ — — ^ gll -)= - r cos 2 <£- r 2 cos 2 T 1 
\ 2 g n o r / 


p 1 _ _1 §_ &24 ~ 1. ^ & 2 4 . _ _ & 24 — „ 772 — P ^ ^ ^ P COS 2 0 

1 _ — 0 n _ ~ 2 f9 " ~~ *> * *2 \ *~ / “ 


24 2g n 5 

1 <9g 




C / 


r 1 = 

33 


2g„ <?r 


13 =_ r (1 + - — £lu)=- r-r^r 1 = L_ T 1 

2 g n d T J 11 cos 22 


r 1 - - 


2*u 5 


d*44 


Yti 1 - 3 /a(7) 2p 2 ( sin ^ + ( i p) 

r + 3 a 

) P 2 (sin 0 ) 

2 m j 
r 

1 - 

7-2 , 2m 

1 + 

T 


2 = 
P 2 (sin 9!)) 





(sin 0) + G~Ty 


r + 3 A(- P 2 (sin^) 1 


4m 

r 


i- 4/ 2 t p 2 ^ sin ^ 


>]) 


and 


12 2 £22 ^ r 


L_ ^22 = 1 + r 


a 


• 2 1 d g 24 _ _ 

14 2 g 2 2 d ** r 2 COS 20 


r 1 

i ^24_ 


r 2 = J— ihz=- tan^+r 1 ; r 2 = ^ ^ 4 =- t , 2 rtan<^r 2 =— 

23 2 g <9<£ 13 34 2 g 22 <90 14 r COS 2 0 


and 


d gn _ _ 1. p 1 . p 3 _ ^ d £33 = _L p 1 = I + r 1 

~ O „ A r 2 13’ 13 9.o rJ r r 2 33 r H 


^ £33 

"gn 
(9 <£ 

1 

^£ 2^2 

^ £33 

(9 0 

1 

d g ?.4 


2 §33 


= — — — — ^ §2 - 2 = — cos 2 0 r 2 = sin 0 cos 0 — cos 2 0 P 1 
a ri rf\ 23 13 


= _ JL_ ^..€24. = _ COS 2 0 r 2 = - tan ii r 1 

24 2 g 33 <90 ^ 34 r r 24 


^Saa -p 1 .2 p 3 • p 3 d _ ^14 = I 


33 2 g 33 <90 *13 *11' ' 44 2 g 33 <90 r 2 34 

and finally 


P 4 = g 24 

12 2*44 


J_ ^^24 _ 1 ^^22 } = _ r p 1 

«24 dT 822 dT 


24 


- - - - - 2 r,\ 

r t 


“ 3 ^4 


44 



r 4 

14 


<? a 


! «44 ^ 




1 ^_g 

2 d ? 


44 


-- r , x i 


r 4 

34 


#24 

2g 44 

1 

2 «44 d( t> 


1 ^24 _ J_ £?? 

«24 ^0 §22 d 4 > _ 


= -^24 


<?g 


44 . 


1 ^^44 = 3*L 

2 d cj> r 


3 h(yJ p 2 ( sin ^ 

— 2 tan — 2 + 2 tan <£j = 0(;n 2 ) = C 

- 4 (- 


^qrV 


sin <£ cos 


These Christoffel symbols will now be substituted into the eqs. of motion (geodesic line) 
_ 4 4 


4!^ 

ds 2 


V £ r • 


Aj 

k = l 1=1 


« S" ^ =° [ z = 1> 2, 3, 4 J 


but first of all the line element ds can be eliminated by the relation 


d*x A Y\ 

ds 2 2j 2j 


j-i 4 /ixA 


Jt= i / = i 

which follows immediately from eq. (223). Therefore 


kl dx 4 dx 4 ys/ 


^ 2 X y (i / 

ds 2 (is \ 


»{1 «o 
1! 

d /dx . 

dx . \dx A 

4 \ 4 

dx 4 \ /dx 4 \ 2 

ds J \ds J 

d 2 x i ix 

ix 2 + dx 

4 

_ ^ 2 *, / 
dx 2 l 

, is 

l>j 2 _ ^ x / 

/ dx 4 

4 4 

22 

k — j 1=1 

r4 dx k 

kl dx A 
4 

dx j 

dx . 

4 

/ ixA 2 
V (is 7 

- 

4 

2 

/c = 

1 2 c 

1 1 = 1 

dx k dx j 

dx . dx A 

4 4 




or after dividing 

by 

(dx 4 /ds) 2 





d 2 X { 

dx} 

4 

■ + 

4 4 

22 

(jl,- 

r 4 dx, \ 
kt dxj 

dxj. 

dx . 

4 

=0 

dx . 

4 


x = ct (x. = c) 


k= l / = i 

Because 

= r ; x 2 = 0 ; x 3 ~ 

the final form of the equations of motion is 

*i + 22 *'=° t«‘- I- 2. 31 

k= 1 i= 1 

The equations of motion in spherical polar coordinates are thus, up to terms of order m, 


r + (I* - 2 r, 4 . ) r 


and 


r,‘, * 1 + r,‘ ■?' . r,‘ r 1 + 2 ( r,‘ s - r,' ) ; ^ ♦ 2 r;, c < 
i * 2 r,V » * 2 r*,» + + 2 (r,V r,‘ 4 f ) . ; ♦ 2 (r,» - r‘.i ) c 4 - 0 


(223) 


(224) 


(225) 


9 = 0 
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<f> + 


r 3 
1 1 1 


22 


e\ 


(r,V- 


2 r 3 4 4 ) 4 > 2 +2uT 3 - r 14 ) r <f> + r 4 3 4 c 2 +2 r 3 4 c e = o 


Using the relations between the different Christoffel symbols and bringing the terms of the unperturbed 
motion to the left side of each equation, the right side will' then be the perturbing acceleration (radial 
component R , lateral component P apd meridional component g, respectively) due to oblateness and 
relativistic effects. Rearraqgemeit of the equations of motion yields 


-^2cos2^-^ 2+ ^ = R m (!!L - r 44 ) C 2+ r,\ 


+ 2 (r; 4 - r,\) r 2 + 2 (r 34 - r^)r<j>-2T^cd 


( 226 ) 


and 


^ ( r 2 # cos 2 0 ) = r 0 cos cos eft — 2 r 6 <f> sin 0 = P 


eft dt 


ee 2 (r 4 4 - r,\ ) r r e cos <f> +2 (r 3 4 4 - r, 1 , ) r Qj> cos +y c C 0 ^ 2 j (r + 2 r cf) tan 0) (227) 


and 


— (r 2 0 ) + r 0 2 sin eft cos eft = r cf + 2r<£-i-r# 2 sin <56 cos cf = Q 


t c r, 1 . 


- - r 4 V + r), j + 2 <r‘, - r,', ) r i i + 2 (r,Vr,‘, ) '4 1 - (r *> W) 


In order to investigate the perturbations of the osculating orbit elements the perturbing ac - 
celeration must be given by its radial components R, its transversal component S , and its orthogo — 
nal component W. Definitions and transformation equations may be found in the author’s report (ref. 
2 5). Introducing y , the inclination of the orbit with respect to the equator of the primary, X. the 
velocity or flight azimuth angle or the angle between the local meridian and the orbit plane, and u = 
co + w, the argument of latitude or the sum of the argument of the pericenter and the true anomaly, there 
exist the following relations for the angles: 


cos <f) sin X - cos y 

cos cf) cos X = sin y cos u 

sin cf> = sin y sin u 


and the following equations for the angle rates 


6 cos 2 cf). = u cos y 

<£ cos cf> = u sin y cos u 


0 cos cf> - u sin X 
cf) = u cos X 


The total velocity is therefore 

v 2 = r 2 + r 2 6 2 cos 2 eft + r 2 <f 2 - r 2 + r 2 u 2 


(229) 


(230) 


(231) 
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Its components 


v r = r * > v 0 = r 8 cos (f) , v ( f ) = r <j> (232) 

are now transformed to the radial, transversal, and orthogonal velocity components, respectively: 


v u = vq sin X + v ^ cos X = r 6 cos <£ sin X + r cos X = r u 
Vy = — v q cos X + sin X = — r 0 cos cos X + r cf> sin X = 0 

The transformation equations for the components of the perturbing acceleration are 
R = R ; S = P sin X + 0 cos X ; W = - P cos X + g sin X 
From the above mentioned relations, there follow the equations 
sin 2 <f> cos X = sin 2 y sin 2 u ; 


(233) 


(234) 


sin 2 sin X = sin 2 y sin u ; 


<^> sin 2 0 = « sin 2 y sin 2 u 


(235) 


which will be used in applying the transformation equations. Substituting the values of P (right side 
of eq. 227) and g (right side of eq. 228) into eqs. (234) and taking the expression tor R (right side of 
eq. 22 6) then (after using eqs. (229), (230), (235) the comp'onents of the perturbing acceleration are 


* = (js - + 2 (r t 4 4 - r 1 1 1 )^ + 2 ^ 2 sin v si n2, 

- 2 (^%)“ cos y 

5 = “ )f 2 sin2 y sin 2 “ + (i]^) T sin2 y sin 2 “ + 2 (r " - r “ )r ; “ 

+ 2 ^ n 4 2 ~| ^- ^ ) T ^in 2 y sin 2 “ + 2 


(236) 


COS20 




(237) 




= - (— j-) — sin 2 y sin u + ( . ^ 3 . — sin 2 y sin u 

\sm 2 0 / r 7 \sin 2 cf> J r 7 

(^) t cos u +- 2 r u sin u ) — ^ 


~ 2 

or, after substituting the values of the Christoffel symbols and writing 

P 2 = P 2 (sin <£) = P 2 (sin y sin u ) ; c - G M = fi 

there is 


(238) 
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IMIIIIII 


- - 
- t t f 1 - 3 '- 



*¥ i' ! [ 1 - : 3 ’‘(rj-] 


and 


5- = 


6 772 y /kV • • . 2 . o 2m r /Q, R\/R\ 

+ / 2 f -yl r a sin 2 y sm 2 a + — 1 f — — Jyyj c u cos 7 


_ Ijl 


2 r 2 
4 772 



in 2 a - (~) 7 sin 2 y sin 2 a 

1 — 3 / 2 r? “ + / 2 ^ r a 2 sin 2 y sin 2 a 


(239) 


(240) 


and 


«' = - | fa h -A(^-) 2 j(f-) 2 sin 2ysin« - | * / a (f)’ 7 sin 2 y sin a 

2772 T -,{Q,R\fR\ C sin y o * * \ 

+ — — r ( )( — J (r cos a +2 t a sm a) 


(241) 


These expressions for the perturbing accelerations can be simplified by using relations for unperturbed 
Keplerian motion in order to obtain first order perturbational effects, namely 


r = — e sm w 

P 


• h 2 (2 1 \ 

; B . -y ; " 2 =^ 7 - -j 


, h 2 [2 \ P 

r 2 = v 2 — r 2 a 2 = v 2 T = « ( — — — — — 0 

r 2 V a r 


. . r (raj a , . » • 2 ( r ^) 2 M { P \ 

r a = 1 = £L («? sin a/) ; r a 2 = = — = A, ±-\ 

r r 2 t r 3 r 2 V/ 

Eqs. (239), (240) and (241) can now be written in the form: 

*- *<•;,(#)* ',♦£{(! - § 

- | / 2 (f/^T ~ T J - 2 7) ?2 “ e sin2y sin sin2 “J 


(242) 


(243) 
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and 


, ft / 3 (R \ 2 . m [ . 

r 2 \ 2 ' 2 V/ r 1 

— J 0 (^) H2 —4 — ^ sin 2 y sin 2zz + 8 e sin w P 2 

2 J 2 \rj \\ s aT / 


and 


W 


/ ft R 2 \ 

3 A / 

ft R 2 \ 2 

— sin 2 y sin 2 J 

r 1 

i\ 

(-) 

e cos y sm zz> + — A y 


1 > 

_ i j 
2 y 2 

l R \ 2 . c, ‘ m 

( — J sm 2 y sin zz + — 

3 

2 

2 \rJ \ a > 

) sin 2 y sin u + 


(244) 


+ 2 T sin y (e sin u/ cos zz +2 — sinzz^+ y A (QJL ^ — sin 2y sin wjj (245) 

The Lagrangian method of the variation of parameters gives the following equations for the time rate 
(variation) of the osculating orbital elements: 


^ ^ [{e sin w)R+(p/r) S] 

at b 


dp _ 2 ft 

z/t h 


t S 


z/e 2 

2e it _ 

1 - e 2 

da 

1 dp 

2? 


( e sin w) 

5T 

2 dt 

a 

dt 

a d t 

6 


ZJ? 77 

dco. 

= — — + cos 
dt 

V fi, “ 2 


P 

(P 


R e sin 

It 

y dt 


h e 2 

LvF 

" V 

da 

dt 

dL 

~~ 71 lit 

dM 

dt 




1 1 
>J bo 


1 

d tt ! 

+ df\ 

71= — 

vr^ 

"e 2 * 

R 


( l -j) 


<7y r cos zz 

zfr h 


W 


daQ r sin u 


dt h sin i 


W 


(246) 

(247) 

(248) 

(249) 

(250) 


Substituting the expressions for the perturbing accelerations R, S, W into the eqs. (246) to (250) that 
yields after reduction: 


da 

dt 


-**(?;(! a© 

m f . . ( 5 3 t\ 3 . (R 
r | \2 4 a) 2 J 2 \r 

■ 3 K 2- «) sin2ysin2 “J -(¥■) 

^— ) sin 2 y sin 2 uj| 


p . 2 

2 e sin zz; . P 2 — — sin y sin 2 u 


(?-«- 


■) 


e sin zz; P„ 


(r)’? ( r+3i (f/ p ?) 


! 4 


(251) 
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1) sin 2 y sin 2 a + 1 4 e sin w 


R' 2 


4 


_ ^ ^2 - ^ - 4 — ^ sin 2 y sin 2 m + 8 e sin w « P J 

f ^CirTf S inV*i»2*l) 


(¥) 


e cos y sin w + 


( 252 ) 


3T - - irr*(t ^(t)1 2 (t - j> p * + (‘ *y ' ,in sin2 ” 

" T [[ 10 (t -i) - (7 — ) (!-«•) f -]- | /.(“/[x (f - l ) ~ 2 (7 +el) 


+ 14 (1 + e 2) L 

P 




6 -(l + e 2 )4 + (l-e 2 ) L 


e sin 2 y sin w sin 2 a 


+ 2 r 


C 1 


OR 2 
h 


(1 + e 2 ) — (1 - e 2 ) T 


r+3 a 


da 

IT 


‘r 

_ 1 A +l) c sin 2 y sin t* sin 2 aj^) 

■vr^ £ . - 2 1 ± {(I - | Z - i)- 


& 


(253) 


(254) 


- ! / 2 (t) 2 [(t- |^- 2 f)p 2 -esi„ 2 ysin I ysin 2 aj 

- 1 <¥)t ^ - K¥T(I) 2 T p [ r+3A (?) 2p |) 

7 T= jh (^T Ja (jJ sin2 y sin 2 “ + 7 |- |/ 2 (t)( 2 - sin 2 y sin 2 a 

+ 2 r (~j~ ) S2n y s ^ n w cos2 u + ■“ s i Q 2 ^ A ^ y sin 2 y sin 2 aj^) 


(255) 


= — (- 3 / 2 -j cos y 


dt pr 


2 u + — ^ - 3/ 2 (^-^(2 - cos y sin 2 a 

+ 2 r ^2? (^2 s * n w s * n 2 “ + 2 j sin 2 sj + 3 A ‘f cos y S2fl2 “Jy* (256) 
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The differencial equations ( 251 ) to (256) have to be integrated with respect to the time or better 
with respect to the true anomaly w using again the relations 

2 

dt — -^r dw , =s 1 + e cos w , u = CO + W 
b t 

of Keplerian motion on the right side of the above-mentioned differential equations. Instead of 
doing this, time mean values of the variable terms over one revolution will be used on the right 
side of eqs. (251) to (256) in order to cancel out short-periodic perturbations. Using 


71 t) 

d M = n dt — — r 2 d w , — ^ 1 + e cos w 

h r 


there is 


f — v sin \ 1 f 

v cos K W ) ~~ 2T 


2 7 T 

-V sin 

t k w dm = 

cos 


r S. 


2 nhp 

The evaluation of this integral for v> k as integers shows that 


n | ,, .v-2 sin 

n _2 f 11 +e cos m;) Kwdw 


=0 (for k >* - 2>0) ; Q^)=0 (forv»2) 


and that 




f 1 for k = 0 
10 k> 0 

1 for k = 0 

1 . 

2 e * = 1 

0 m > I 


e 2 

1 + 2 “ for k = 0 

e k - 1 

1*2 K= 2 

4 

0 K> 2 


1 + e 2 for k = 0 



0 k > 3 
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1 + 3 e 2 + J- e 4 for k = 0 
o 



L 0 k > 5 

Due to u = a) + the addition theorems 


sin ^ • sin 2 z/ = — ^ cos (2 co + 3 u>) + cos (2 co +u>) 
sin w * cos 2 « = — sin (2 co +3 w) — sin (2 co + w) 


and the expression for the second harmonics 

_ i I7i_ i 

2 Lv 2 


= -j- (3 sin 2 


sin z z/ 


i) 




cos 


52 




there is 




sin w sin 2 u ' 


sin tu cos 2 u 


f sin w cc 

V ^ 



^ 0 for v > 3 


These general formulas provide now the following special time mean values applied to eqs. (249) to 
(254): 
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/sin w sin 2 u\ 3 n /, e 2 \ 0 

( Fi > 4 Fpa e (? + 6 >° s2w 


Substituting these time mean values into eqs. (251) to (256) then, after reduction, the equations can 
be written 


a 45 (m \ /R V l+e?/6 

V T "°\Fj ’ W 


, 2 0 


e, 2 sin 2 y 0 sin 2co 


p 3 /W2\_/^\ 2 o.9 . 0 

-o M~ r )j 2 ( — ) ^ 0 2 sin 2 y 0 sin2 ^ 
Po 2 \ Poy \ PqJ 


(257) 


(258) 


p \ 3 (m \ ( R \ 2 2 A o ^ 2\ , 2 . c 

_ = 2e 0 e= (l-. 0 no (-Jj 2 (-Je 0 (l^ +J e 0 Js ln y 0 S1 n2m 


thus 


* = " 5 % ^ 3+ ^ e o^) si " 2 yo s i“2a, 


* = “ + a *f2 cos > / = ”o(f ^07)(l" 1 Sin2 ' 


- - / 2 
2 72 w 


(259) 



■ y oj + 


3 - 


16 + ?j- e^)(l - | sinVo)+ ^-(13+23 eg) sin 2 y 0 cos 2 a, 
. p /qr 2 \ 3 A Air 2 Y/. 3 . , n 


(260) 


1 ■ "4 KfX 2 ■ f f KfXXX'X - f sin Y 

— ( 1 + 2 e Q 2 ) (2 - | sin 2 y^) + i (13 + 23 e 2 ) sin 2 y Q sin 2a> +-i-e Q 2 cos 2 y Q cos 2a, 


- 6 r 



r o 



(261) 


I /2 (t - ) 2 t ( x + t e °) ( 7 ~ 4 sin2y o) “ i (13 + 37 e o ) sin2 y 0 cos 2 ® 
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y =- 3 « 


3 (m \ , (R V „ . „ 

° (ir ) M t e o Sln 2 y o sin 2 1 

\"o/ \ p 0/ 


3 . /R \ 2 


to \ 13 . / R 


a “ ”°\2" ; 2 (yj cos *> + ("pJ|2 M pj C0S *> 1 + 2 e o~ 2 e 0 cos 2 


or/OR J \ 3 A /flR 2 ^ 2 

- 2r ljrr 2 A hJ '°‘r 


It is easy to integrate these equations where the orbital elements on the right side have their unper- 
turbed constant values (subscript 0) except co, which is a linear function of time t, namely 


*”o! '=( 0 ( 2 - | 


(i) n + CO • , ) 

0 non-rel. v 0 ' 


The integrals 


I COS 2 CO dt = ~r 


Isin 2 co dt =s — 


show that long-periodic perturbations appear in addition to the secular perturbat 
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B. The Time-Dilatation Effect 


Using the line element of eq. (222) the ratio of the proper time element d T to the local time ele- 
ment dt is given by 


dr\ 2 / d s \ 2 _ , _ 2 7?? f, _ /Re 

dt 7 ~VdUj ~ ~ \ J >\rs 

-{ 


2 , . ,, (ClRe' 2 

P (sin <p) +( 


r+A (~) P 2 (sin <f>] 


+ 0 (m 2) - |l + 1 - ] 2 (^j P 2 (sin 0) 


1 v 2 2 m 2 , dQ 

JT 2 + T~ r \c^ / cos ^ Tt 


or accurately enough 


tfr\ 2 = 1 _ 2w 


dt 


r L 


1 + | / 2 (t^ (1 “ 3 sin2< ^ 


thus 


d r 
dt 


i- ^ 


i + 


i ; 2 (fJ(l-3sin 2 0) 


v 2 t; 2 /2 + U 

7T~ — -L O 


(265) 


This formula, containing the invariant proper time element dr — ds/ c, is now applied to an earth clock 
(r = R, v = Q R cos 0) : 


dr „ l _ HL 
dT E - " R 


1 + £ 1 2 (%-) d - 3 Sin2 ^ 


l n 2 /? 2 cos 2 </> n 

— - « “ ~ 1 — o “ 


const. 


Assuming the Earth’s surface an equipotential surface, the potential W Q is constant and can be expressed 
by its value at the equator, thus 


1 

1 

/QRV , 

m 

\ i . i 

1 + 2 ^ 2 

2" 

( o) - 1 - 
V c y 

K 

[ 1 + 2 ^ + 2 X 


dr , 77i 

TtS " IT 


where y = Q 2 R 3 /\i has been introduced. 

Applying eq. (265) to a satellite clock moving with the satellite in an arbitrary orbit 
(sin = sin 8 = sin y Q sin u ; v 2 /c 2 = m ( 2/r — l/^ 0 ) > u = co +w) ? 

that yields 


( 266 ) 


r i 

f R \ 

y + 2 

(v 


, 1 m m 

l + T> 

2 r 


(1-3 sin 2 y sin 2 uy 

1 / d \ 2 / 1 

2 + ± / 2 (^) ^1 — 3 sin 2 y 0 — 


m (2 _ 1 

2 ir a 


2 « 


Using the time mean values 


1 ; (l°j 3 = (1- e 2 ) 3/2 ; -^cos 2 
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then results 



i; H Gt)’' 1 -'”'" C 1 - * * , 4 


(267) 


Dividing eq. (266) by eq. (267) gives for the relative difference in the time rates of a satellite clock 
compared wi th a standard earth clock 





+ 




(263) 


Eq. (268) consists of three main terms. The first term (gravitational red shift and time dilatation) was 
first given by Winterberg (1955, Ref. 26) and Singer (1956, Ref. 27). The former author also added the 
last term due to the rotation of the Earth. The second term, due to the oblateness of the Earth, corrects 
and generalizes the term given by Hoffmann (1957, Ref. 28) for circular orbits in the equatorial plane.. 


This paper will be concluded by listing certain constants of the Earth’s gravity field which appear 
in the text, namely 


*. 

h 

r 

A 

m 

m/R e 

n 

a R e 
SIR* 

x 

a R e /c 
(Q,R e /c ) 2 


light velocity in vacuo, 299792.50 + 0.10 km/sec (K.D. Froome, 1958) 

gravity factor of the Earth, GM = 398613-52 km J /sec 2 (Herrick, Baker, 1957) 

equatorial radius of the Earth, 6378.150 + 0.050 km (“Baker, 196 1) 

oblateness constant of the Earth, (1082.190 + 0.023) x 10" 6 (Kozai, I960) 

inhomogeneity factor of the Earth, 0.3336 ~ 1/3 

inhomogeneity factor of the Earth, — 0.0429 ~ — 3/70 

gravitational radius of the Earth, fi/c 2 = 0.4435 cm 

potential energy factor of the Earth, 6.95377 x 10 ~ 10 

angular velocity of the Earth’s rotation, 7.292115083 x 10 -5 sec- 1 

rotational velocity at the Earth’s equator 465.102 m/sec 

angular momentum (per mass unit) at the Earth’s equator, 2.96649 x 1 0 9 m 2 /sec 
centrifigal factor at the Earth’s equator, O 2 R$ /\i = 3461.30 x 10“ 6 
rotational velocity ratio, 1.55142 x 10' 6 
square of the rotational velocity ratio, 2.40691 x 10 “ 12 
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